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We study the possibihty that particle production during inflation could source observable gravity 
waves on scales relevant for Cosmic Microwave Background experiments. A crucial constraint on 
such scenarios arises because particle production can also source inflaton perturbations, and might 
ruin the usual predictions for a nearly scale invariant spectrum of nearly Gaussian curvature fluc- 
tuations. To minimize this effect, we consider two models of particle production in a sector that is 
only gravitationally coupled to the inflaton. For a single instantaneous burst of massive particle pro- 
duction, we flnd that localized features in the scalar spectrum and bispectrum might be observable, 
but gravitational wave signatures are unlikely to be detectable (due to the suppressed quadrupole 
moment of non-relativistic quanta) without invoking some additional effects. We also consider a 
model with a rolling pseudoscalar that leads to a continuous production of relativistic gauge field 
fluctuations during inflation. Here we flnd that gravitational waves from particle production can 
actually exceed the usual inflationary vacuum fluctuations in a regime where non-Gaussianity is 
consistent with observational limits. In this model observable B-mode polarization can be obtained 
for any choice of inflaton potential, and the amplitude of the signal is not necessarily correlated 
with the scale of inflation. 



Contents 

I. Introduction [2] 

II. Formal expressions for observable correlators H] 

A. Decompositions, and sourced equations for h\ [5] 

B. Formal solutions and correlators [S] 

C. Phenomenology [7] 

III. Model I: Vector produced by non-adiabatic change of its mass [5] 

A. Scalar perturbations sourced by the vector modes [TO] 

1. Master equation for ( [TO] 

2. Scalar source and power spectrum [11] 

3. Scalar bispectrum [TO] 

B. Gravity waves sourced by the vector modes [13] 

C. Phenomenology [TO] 

D. Comparison with GW sourced by modes of different spins [TO] 

IV. Model II: Vector produced by a pseudo-scalar interaction [TO] 

A. Scalar perturbations sourced by the vector modes [21] 

1. The master equation [5T] 

2. Two-point and three-point correlation functions [22] 

B. Gravity waves sourced by the vector modes [22] 

C. Phenomenology [23] 

V. Conclusions [26l 

Acknowledgments [25] 

A. Longitudinal vector mode [25] 

B. Source for ( in Model I [3l] 



2 



C. Source for in Model II 



D. Fermionic production and gravity waves 



m 



References 



I. INTRODUCTION 



Inflation is currently the standard paradigm for solving the flatness, horizon and the unwanted topological relic 
problems in standard big bang cosmology. In addition, the quantum fluctuation of the inflaton field provides the seeds 
of inhomogeneity that seem to fit very well with the current observation on the Cosmic Microwave Background (CMB) 
and Large Scale Structure (LSS). A particularly important probe of the physics of infiation is provided by primordial 
gravitational waves (tensor fluctuations) . A number of observational searches for gravitational wave perturbations are 
proposed or are currently underway; these may be probed through the B-mode polarization of the CMB [IHS] or, on 
smaller scales, through interferometers such as LIGO H], VIRGO [S], DECIGO jS], Einstein Telescope [7], or LISA 



During inflation, gravitational wave fluctuations are inevitably generated by quantum fluctuations of the tensor part 
of the metric. These have an amplitude controlled by jjj, where H is the Hubble scale and Mp « 2.4 • 10^^ GeV. The 

gravitational wave signal from vacuum fluctuations is detectable only when the inflaton fleld range is trans-Planckian 
[S], which might be challenging to realize in a controlled effective fleld theory. 

Additional sources of Gravitational Waves (GW), which are uncorrelated with the usual quantum vacuum fluctu- 
ations, may be present in the early universe (see, for instance, |10j for a recent review). Two broad categories of 
mechanisms are: 

1. Models involving phase transitions: For example, in first order phase transitions, vacuum bubble collisions 
pn [T^ and the subsequent turbulence [13] could source GW. The generation and decay of cosmic strings can 
also give rise to large GW 14J. Along similar lines, the self-ordering of a scalar field after a second order phase 
transition has also been considered 15 . 

2. Models involving particle production: Generically the inflaton should be expected to couple to some 
additional degrees of freedom, as would seem to be necessary for successful reheating P^20j . In this case, 
there is a natural possibility that the time-dependence of the inflaton condensate during inflation leads to the 
production of some other degrees of freedom which may, in turn, provide an important new source of GW 
pm23j . A variety of different models have been proposed; see below for more discussion.^ 

Given the importance of gravitational waves as a probe of inflation, it is important to understand if such mechanisms 
could be competitive with the usual spectrum of GW from vacuum fluctuations. ^ 

Any mechanism which is being invoked to source gravitational waves might also source scalar metric perturbations. 
Therefore, one must take care not to spoil the usual prediction of a nearly scale invariant spectrum of Gaussian scalar 
curvature perturbations. Sometimes this concern is evaded by restricting attention to effects that take place on the 
small scales relevant for interferometers where the scalar fluctuations are not strongly constrained. In this work we 
will mostly be interested in models of particle production during inflation, where gravitational waves are sourced on 
CMB scales. To test the feasibility of such scenarios, it is crucial to study also the spectrum and bispectrum of the 
scalar fluctuations, to ensure that these are consistent with observations. 

Models of particle production during inflation have received considerable attention in the literature; see for example 
[211f23l I34ti48] . One class of models involves instantaneous bursts of particle production, leading to localized features 
in the cosmological perturbations. Early studies focused on the production of fermion [37] or scalar [3H] particles, 
neglecting the feed-back of the produced quanta on the perturbations of the inflaton. In Ref. [35^ it was however shown 
that this feed-back effect actually dominates observables for the case of scalar particle production. Ref. 39 considered 
a simple model where scalar x particles are produced through a coupling g^'^P'x^ to the inflaton, using both lattice 
fleld theory simulations and also analytical methods. Models of this same type were subsequently analyzed in the 



^ For gravitational waves from particle production at the end of inflation, see I24H28I : effects on the scalar fluctuations were discussed in 
[29H32]. 

^ A GW signal may also be left at the largest scales as an imprint of a pre-inflationary era if inflation had only a minimal duration | 33| . 
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context of trapped inflation f?7^ . More recently, models of instantaneous vector particle production have been studied 
in connection with GW at interferometer scales [21j . Several scenarios of GW from scalar and string production were 
instead studied in f^ . 

Another possibility is that particle production occurs continuously during inflation. This is quite natural in the 
context of axion inflation [44, . Even in the simplest models of inflation driven by a single axion in slow roll on a smooth 
flat potential, pseudoscalar couplings LpFF to gauge fields are ubiquitous. This interaction leads to a continuous 
tachyonic production of gauge field fluctuations during inflation. There are a host of interesting phenomenological 
signatures: observable equilateral non-Gaussianity [111133115], GW at interferometer scales t21j|23j, and excess power 
at small scales [501 EI]- Additionally, backreaction effects in such models can assist inflation by dissipating the 
kinetic energy of the inflaton I^Sl W3 , analogously to warm inflation '34] and trapped inflation |47] . ^ In [54lf56] a 
model of gauge field production through the scalar coupling l'^{ip)F'^ was considered, in connection with primordial 
non-Gaussianity and, perhaps, magnetogenesis (this last application is problematic |56[ 157) ). 

Above we discussed some simple models of particle production during inflation where quanta are produced by 
a direct coupling between the inflaton and some additional degrees of freedom. In such a scenario the produced 
particles interact with the inflaton through couplings that are typically much stronger than gravitational. Hence, 
these produced particles will tend to source scalar curvature fluctuations much more efficiently than GW. It is not 
surprising, therefore, that the most stringent CMB constraints on such models often come from features and non- 
Gaussianity in the scalar fluctuations, rather than from GW. 

In this work, we investigate the possibility to source a significant GW signal on CMB scales without ruining the 
spectrum of curvature fluctuations. To minimize the impact on the inflaton fluctuations, we assume that particle 
production takes place in a "hidden" sector that is only gravitationally coupled to the inflaton.'* We compute for 
the first time the scalar perturbations induced by the gravitational coupling, and we compare their phenomenological 
impact with that of the produced GW. Since gravitational interactions are unavoidable, the case we investigate can 
be thought of as a kind of "best case" scenario for the production of GW while minimizing the effect on the scalar 
fluctuations. Concretely, we will focus on two models involving the production of spin-1 particles: 

1. Model I: Particle production takes place in a hidden sector where a local U{1) invariance is spontaneously broken 
by the expectation value of some complex scalar ip. Here gauge flelds are produced at an isolated moment - 
namely, when the mass of the gauge fleld crosses zero during the evolution of ■0 - and we have localized features 
in the scalar spectrum and bispectrum, in addition to a localized feature in the tensor spectrum. For a single 
isolated burst of production, we find that observational constraints on the scalar perturbations exclude any 
interesting effect in GW. However, in a concrete model there may be many bursts of particle production and 
their resonance could enhance the effect. We note that this model could lead to interesting phenomenological 
signatures in the scalar fluctuations; localized features in the spectrum and bispectrum are discussed. 

2. Model II: Particle production takes place in a hidden sector where a rolling pseudoscalar sources gauge fleld 
fluctuations continuously. Here we flnd that GW from particle production can be competitive with the vac- 
uum fluctuations, or even larger than the vacuum fluctuations, without violating observational bounds on non- 
Gaussianity. In this model the primordial tensor spectrum can be detectable for any choice of inflaton potential 
(that is, also for small field inflation). In this model, GW from particle production are chiral. We show that 
parity violation in the tensor sector can be almost maximal, which provides a distinctive observable signature 
of this scenario; see also [^ [5^ . 

In summary: we find that the possibility to source interesting GW from particle production is rather model 
dependent. Even if particle production occurs in a hidden sector, coupled only gravitationally to the inflationary one, 
this is not a sufficient condition to ensure that observable GW can be sourced without ruining the scalar spectrum. 
On the other hand, we find that in some models GW from particle production can actually exceed the usual vacuum 
fluctuations. The reason for the smaller GW production in Model I with respect to Model II is that in Model I 
the gauge quanta are highly non-relativistic after their production, which suppresses their quadrupole moment. We 
verifled that the same suppression takes place if, instead of gauge flelds, the particles produced in this mechanism have 



^ The pseudo-scalar interaction is also used to [53 to dissipate the inflaton kinetic energy from its classical interaction to a non-abelian 
vector field. 

* Our scenario differs from the curvaton model [58] . in which the inflaton provide the energy density to support inflation while a light 
curvaton fleld provides the seed of scalar perturbation. The GW spectrum in the curvaton model has been considered at tree-level 
1591 1601 and at 1-loop level |61| , with possible detectability at future interferometer experiments |61| . See |62| for examples of vector 
curvaton. 
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spins or i. ^ As we mentioned, a greater GW effect may be obtained for multiple instances of particle production, 
or for a more complicated evolution of the gauge field mass, or if the massive particles decay into massless ones short 
after they are produced P^ . 

Therefore, due to the possibility of GW from particle production, a measurement of primordial B-modes does not 
necessarily constitute a measurement of the scale of inflation. Nor does a detectable B-mode signal necessarily require 
super-Planckian excursions in field space. Fortunately, as we discuss below, the produced GW may be distinguished 
from the vacuum signal, either by their localization at some given wavelength (assuming that Model I can be modified 
so to enhance the GW signal) , or by their violation of parity. 

In terms of observational prospects of gravitational waves on CMB scales. Model II is certainly more promising, as 
the gauge quanta are automatically relativistic and their quadrupole moment is not suppressed. While the focus of the 
present work is to examine the signatures and constraints of gravitational waves produced during inflation as a result 
of hidden sector particle production, the UV sensitivity of inflation motivates us to ask whether such models can be 
realized in string theory. Indeed, the 4D low energy spectrum of string theory contains a myriad of axion-like fields, 
e.g., those arise from the reduction of antisymmetric form fields on cycles of the internal space. These closed string 
axions couple to U{1) gauge fields on the worldvolume of D-branes via ipFF couplings. The axion decay constant 
/ one typically finds in string theory constructions is of the order of the GUT scale Mqut ~ 10^® GeV (see, e.g.. 



[63ll64) ) which sits comfortably within the allowed window for consistency of our model (see eqs. ( [115^ and (116)). It 
is not difficult to arrange the inflaton to have no direct coupling to the axion and the gauge field. We outline some 
ideas to realize Model II in string theory in the concluding section. We leave, however, a detailed study of these and 
other string theory embeddings for future work. 

This paper is organized as follows. In Section [ll] we develop a very general formalism which can be used to compute 
the effects of particle production on the scalar and tensor cosmological perturbation in a variety of models. In Section 
|III| we apply this formalism to Model I, discussed above. In Section |IV| we instead apply our formalism to Model 
II, discussed above. In Section [V] we conclude. In Appendix |Xj we discuss the production of the longitudinal mode 
in Model I. This is the first time the longitudinal mode has been accounted for in a model of production of massive 
vector fields during inflation. In Appendix [Bj we discuss some details of the computation of the scalar perturbations 
for Model I. In Appendix [C] we perform an analogous computation for Model II. In Appendix [Dj we discuss the 
production of massive fermion fields during inflation, and the amount of GW that they generate. 

II. FORMAL EXPRESSIONS FOR OBSERVABLE CORRELATORS 

We assume that quanta of a vector fleld are produced (either at some discrete moment , or throughout inflation) 
in a sector only gravitationally coupled to the inflaton. We are interested in computing scalar and tensor cosmological 
perturbations sourced by the vector quanta. The vector field can be decomposed as 

Ao = , A, = y -^e^*^ - (r, fc) , (r, fc) = ^ ej^^ (fc) [ax (fc) Ax (k) + 4 (-fc) A^ (fc)] (1) 



A=± 



where we have included only the transverse vector polarization (in Section III we actually study also the longitudinal 
polarization, as it is present, and may be relevant, in that model). Here are circular polarization vectors satisfying 

k ■ e**^*) (jtj — 0, k X e**^^^ (^kj = ^ike^^^ ^fc^, e (~^) ~ (^) ' ^^'^ normalized according to 
eW (fc)*.r(^') {k)=Sxx'. 

In this section, we introduce the formalism for the correlators of the scalar and tensor perturbations sourced by 
this vector vector, under the assumption that the latter is only gravitationally coupled to the infiationary sector, and 
that it enters quadratically in the energy-momentum tensor. The discussion is divided in three parts; in Subsection 



II A we present the structure of the equations for the scalar and tensor perturbations in presence of the source; in 
Subsection |II B| we present the formal solutions, and the structures of the correlators; in Subsection |II C| we show how 
these expressions enter in observable quantities. 



^ We find that the resuhs agree in all three cases, with only order unity differences coming from counting the number of degrees of freedom 
and also from spin statistics. These order unity effects can be encoded in a very simple formula, which we present. 
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A. Decompositions, and sourced equations for 



We compute scalar cosmological perturbations in the spatially flat gauge (Sgij. scalar = 0. In this gauge, the curvature 
perturbation is (up to negligible corrections, as we will see in the next two sections) C = S(p, where and 

S(p are the unperturbed and perturbed part of the inflaton, respectively. ^ We decompose 



(2) 



As we show in the next sections, the effect of the vector fields on the inflaton modes can be encoded in the 
approximate equation 



where the source is formally of the type 

Jf (t, fc) = 



3/2 ^V.i] 



(j, kj ~ (j, kj 



(t, k, jjj A, (r, p) A-j {t, k~fj 



(3) 



(4) 



where O is a model-dependent operator. 

In the tensor sector, the gravity waves are encoded in gij — a} (Sij + hij), where the modes hij are transverse and 
traceless. We introduce the canonical modes 



Mpa, 



which obey the equation 



d^k 



Ak 



(2^) 



3/2 



^ ^ n,,,, (fc) {k) , n,,,, (fc) = ' (fc) ef^ ( 



(5) 



(6) 



The source is obtained by starting from the transverse and traceless spatial part of the energy momentum tensor 
and projecting along the A polarization with One finds 



/d^ X 



-ik-x 



^T,,(r,f) 



(7) 



We note that, the multiplication by 11*^ ^ automatically projects on the transverse and traceless part of Tij. If (as in 
the cases that we will consider) the energy-momentum tensor is quadratic in the gauge fields, we recover an expression 

formally identical to (4), with a different (and model-dependent) operator Ox^ij [t, fc, pj . 
In summary, we have formally identical equations for the scalar and tensor canonical modes 



a 
a 



(r, fc) ~ Jx (r, fc) 



(27r 



,3/2 



O 



X.ij 



(t, k, pj A, (t, p) Aj (T,k-pj , X^{ip,X = +,X = -} 



(8) 



The operators Ox depend on what X is, and also on the model. In all cases, they are invariant under the simultaneous 
« -H- j and p^k — p operations. 



Throughout this work we denote physical time by t and conformal time by r. Derivative wrt physical (conformal) time are denoted by a 
dot (prime). We denote the scale factor and the Hubble rate by a and H, respectively. We denote by Mp the reduced Planck mass. We 
use the mostly positive —,+,-!-, -|- signature everywhere in the paper apart from the model of fermionic production (Subsection |III P] 
and Appendix |D| where we switch to the opposite signature. 
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B. Formal solutions and correlators 



The equation ([s]) is formally solved by 

Qx (r, fc) = Qx,v (t, fc) + Qx.s (r, tj (9) 
where Qx.v is the standard vacuum solution, obtained in absence of the vector source; the sourced term is instead 

(r, fc) = r dT'Gu{T,T') Jx{T',k) 



Gk {r, r') 



— [kr cos (fcr') — sin (fcr')] 



-fcr < 1 



(10) 



where the Green function is only given in the super-horizon regime —fcr <C 1. The two terms in ^ are uncorrelated. 

We are interested in the two and three point correlators of the sourced solutions. For the two point correlator, we 
find 



where 



Qx,s (t, Qy,s (t, ^2)) ^ J driGki (r, ti) J dT2Gk^ (r, T2) (jx (n, ^i) Jy {j2, h 



(11) 



X I A 



(n, pi) Am 



T2, k2 -P2 



) (i, (n, fci-pi)i,(^2,P2)) (12) 



In writing this expression, we Wick decomposed the (A*) correlator coming from ( J^) , and then used the symmetry 
of the Oyim operator. 

Prom Q we then see that the two point correlator of the gauge field is formally of the type 

(13) 



(n, ql) A, (t2, 92)) = Vlp (gi) 2?;o%) [n, T2; gi] 5(3) (gi + q2) 



where 



(14) 



The index (0, 0) on 2? indicates that no time derivative is acting on the modes A on the l.h.s. of ( 13 1, and it is introduced 



for later convenience. Using the hermiticity of (x), one can show that f |q p) ['''ii ''2; l] =^ -^(o 0) [^2, ti; q]- We show 
in the following sections that, in both models we study, 2?|o 0) bscomes real and symmetric under ti -s-)- T2 for all times 
after the particle production. 

Inserting (131 into (12 1 we obtain 



Jx {n, h) Jy {r2, k2)) = 2<5(3)(fc, + fc,) J £P ^^^tJ {p)vlp (f^^i) 

Ox,t3 (ti, ^1, p) (t2, -fci, p - fci) P^ol,) [ti, T2; p] 2?[o Ti, T2; |p- fell 



(15) 

Finally, we only need to compute the connected three point correlator of Q^p. Proceeding as for the two point 
function we obtain 



nQv'.s (j, ^*)^ = y" '^■^i^fei ('^' "^i) J dT2Gk2 (t, T2) J drsGks (r, T3) ^11 (j^^ ^i) 



(16) 
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with 



J if ( ''"i I 



\i=l 



(2^0 



9/2 



(^Ti, fci, O^^lm (t2, fc2, P+ ^2) O^.no (^3, p-ki^ 



^(0,0) I'^i' '^2; p] 25[o.m '^1' '^3; |p - ki 



(0,0) 



(0,0) 



T2, T3; Ip + fczl 



(17) 



C. Phenomenology 

From the two point scalar correlator we obtain the power spectrum 



(18) 



The two solutions (|9| are incoherent and therefore their power spectra add up. Recalling the standard slow roll result 
for the vacuum mode, and using the two expressions (11) and (15 1, we obtain 





= PC.V+PC.S 






PC,V 


47^2(^(0)2 




kl i/2 


Pc,^ 


7r2 a20(")2 



fdnGk, (r, n) fdr.Gk, (r, r,) J ^ (P^T'iP {p^^) 

J {t1, kl, pj d^,lrn ( 



i2nr 

(0,0) 



T2, -kl, p- kl ) 2?!^],^ [ti, T2; p] "^[q^) n, ^2; \p- ki\ 



(19) 



For the tensor mode, we show below that modes of different polarizations are uncorrelated. Starting from the 
decomposition ([5|, the power in each polarization is 



kf 



27r2 a2M2 



Qa (kl) Qx (fc2 ) ) = Px S^'^ (kl + k2 



(20) 



Also in this case, the power of the vacuum and the sourced modes add up 
Px 



Px,v 
Px,s 



Px.v + Px,s 
i/2 



7r2 Afp2 



(27 

OA,y (ti, kl, pj dxdm (t2, -kl, p- k^ '^%]o) [n, T2; p] P^gJ) Tl, T2; |p- ki\ 



The tensor-to-scalar ratio r is defined as 



(21) 



(22) 



If the sourced term is absent, one recovers the standard vacuum slow roll result Tv = ^ j^j2 — 16e. However, we see 
that the sourced contribution can modify this result. 

Finally, we are interested in the bi-spectum of the scalar modes, defined as 



C (fci) C (^2) C (4) ) - 5c {k^) ^''^ {kl 



(23) 
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We disregard the contribution from the vacuum modes, as it is unobservable for standard slow roll inflation. Using 
(161 and (17) we obtain 



B^iki) =8 driGk^T, Ti) / dT2Gk^{T, T2) / dT^Cksir, T3) 



d^p 
(2^) 



9/2 



(7,a' ^a" 



(0,0) 



T-2, T-3; |p + fc2 
(24) 



In the following sections, we compute the bispectrum for equilateral configurations (which is where it is peaked in 
the models that we consider), and use it to define an equilateral nonlinear parameter 45J 



/7VL,cquil. cff. (fc) 



10 l_ R I 

s5/2 p2 fj^\ ^C\ki=k2=k3=k 



9(2^)5/2 p2(fc) 



(25) 



III. MODEL I: VECTOR PRODUCED BY NON-ADIABATIC CHANGE OF ITS MASS 

We consider a model where the sector in which particle production takes place has a local U(l) invariance sponta- 
neously broken by the expectation value of a complex scalar '5: 



S — d X\J 



1 1 

-fR - -(d^f - V{^) - I (5^ - ^eA^) M^p - C/ (|*|) - -F^ 

" V ' " V ' 

infiaton sector hidden sector 



(26) 



The field (/? is the infiaton field, which is assumed to be only gravitationally coupled to the — sector. We work 
in the unitary gauge ^ = where ■(/; is real. We assume that the background value V'^'^-' (i) crosses zero at the time 
t* during inflation. Close to this time, the gauge field mass can be approximated by 

m = e ?/'^°' e^'f'^ {t - U) = rh* {t - U) (27) 



For definiteness, we take m* > 0. The (comoving) frequency of a gauge mode, oj — \/k^~^~a?rnP varies nonadiabatically 
(w' > w^) when the mass vanishes. As we will show later, the gauge field modes that dominate the observational 
signatures have k ~ Vrh^. For such modes, the frequency changes nonadiabatically in the time interval ^ t^zL yi=, 

provided that the expansion of the universe can be disregarded in this time interval. This is the case for ^ -g.. In 
the remainder of this Subsection we concentrate on the production during this interval, and disregard the expansion 
of the universe (we use physical and conformal time interchangeably). We can eliminate any reference to the scale 
factor by normalizing a (i*) = 1. In the following Subsections the expansion is taken into account. 

The non-adiabatic change of the frequency causes non-perturbative production of the gauge modes. To compute 
this, we decompose the vector field as in ([T]). We further decompose the mode functions in positive and negative 
frequency modes 

A^{k)^ak{T)fk{T)+Pkir)rk{T) , fk^ "" '^J^ - 

^2w (r) 

A'^ (k) = -iLO K (r) fk (r) - Pk (r) fk (r)] (28) 

(the second line is the decomposition of the modes of the conjugate momentum to Ai). The decomposition ([T]) 
disregards the longitudinal vector mode; we compute this mode in Appendix and we discuss its effects below in 
this section. We have suppressed the index A in this decomposition, since the Bogolyubov coefficients are the same 
for both helicities. For t <^t^,, the mode is in the adiabatic vacuum a = 1 (up to an arbitrary phase); the quantity 
|/3p is the occupation number of the gauge modes. 



The gauge field modes satisfy A" + uP'A = 0. With the approximated expression (27), this equation can be 
analytically solved in terms of two parabolic cylinder functions. The linear combination satisfying the proper initial 
condition a = 1 gives (up to an arbitrary phase) |171 121j 



afe(t>i,) ~ Vl + «^^ , /3fc(i>i*)--e""^ (29) 
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Namely, one finds a O (1) occupation number for modes up to fc ~ \fm*i while the production is exponentially 
suppressed at higher momenta. An analogous result is obtained for scalar [17j and fermion (37j fields produced by a 
mass varying as in (pTl. This result is valid provided that 



(30) 



where (as we mentioned previously) the first condition ensures that the expansion of the universe can be disregarded 
during the interval of particle production, while the second condition imposes that the kinetic energy of tp is negligible 
with respect to the infiaton energy density. 

Having the mode functions, we can now compute the correlators (13). The correlators need to be regularized. 
The regularization can be performed by normal ordering with respect to the time dependent annihilation creation 
operators 



a\{T, k 



ak (t) ax (fc) + PI (t) al (-fc) 



l3kiT)ax(k) +al (r)4 



(31) 



These are the "physical" operators of the system, as one can show that they diagonalize the Hamiltonian at all times. 
Using ([T]) and (28) it is immediate to show that, in terms of these operators. 



(r, fc) = E ''i'^ (^) [fk (^) ^) + A* (^) «1 



(32) 



We then compute ^: Ai ^r, kj Ai [t' , k'^ :^ by normal ordering with respect to the a, and a\ but by recalling that 

the vacuum state is annihilated by the original time-independent a (we are working in the Heisenberg picture, in 
which the states are constant). ^ Proceeding in this way, and casting the result as in (13 1, we obtain 



^[o!o) k] = ak (r) (r') /, (r) /, (r') + (r) a*, (r') /,* (r) /,* (r') 

(r) PI (r') A* (r) /, (r') + PI (r) P, (r') f, (r) A* (r') 



(33) 



Clearly, the result is cj— independent, since both helicities are produced in the same amount. After the particle 
production, this expression simplifies into 



V\^l^ [r, r'; k] = a, PI h (r) h [r') + «I. Pk ft [r) ft W) + [/* (r) /, (r') + /, (r) /* (r')] 



(34) 



where a and P assume the asymptotic values ( 29 ) . As we anticipated, this expression is real and symmetric under 

without normal ordering, and then sub- 



T O r . 

The same result (34 1 is obtained if one computes I^Ai 

tracts the term that one would have in absence of particle production (a = 1, /3 = 0), as done in [39]. Therefore these 
two regularizations are equivalent after the particle production has taken place. 
For further convenience, we extend the definition ( 13 ) to 



d_ 



At (n, qi) 



d_ 



)' A, (r2, q,) :\ ^ ^ (^O^i) [n, r,- q,] 6^'^ (q, + q^) 



(35) 



We compute these quantities by using (28). After a and P have assumed the asymptotic values (29), the correlators 
of our interest become 



V}';;,^ [r, r'; k] = -zu^ (r) {a, PI h (r) h (r') - al Pu ft (r) /,* (r') + \Pu? [- fk (r) fk (r') + A (r) /,* (r')]} 
2?[o,\) [r, r'; k] = -ic, (r') {a,. Pl h (r) h (r') - P^ f* (r) /* (r') + [/* (r) A (r') - A (r) A (r')] } 

[r, r'; k] = -u^, (r) o.^ (r') {a^ /3;5 A (r) A (r') + /3, /,* (r) /,* (r') - \P,\^ (r) A (r') + A (r) /,! (r')] } 

(36) 



This procedure is conventionally adopted in Bogolyubov computations, although it is not always spelled out; indeed, the initial Hamilto- 
nian, without normal ordering, can be cast in the form H = f dP'kuJi^ N, where the counting operator is N = 2£__t£_a One then defines 
normal ordering wrt the barred operators, ; iV := a, and evaluates (: N :) using \3l\ and recalling that the vacuum is annihilated by 
the original time-independent a; only in this way one obtains {: N :) = |/3p. We employ the exact same procedure for evaluating the 
correlators D. 
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We note that in the regime in which ( 34 ) and ( 36 ) are vahd, Ot'd'"^ 



(0 0) ['''' ''''' ^] ~ "^(10) '^l' ^^'^ analogously for 



the other terms. Namely, in this regime, time derivatives can be equivalently taken before or after evaluating the 
correlator. One can verify that (for the regularized correlators) this is not the case during the particle production, 
when a and /3 are still functions of time. 



A. Scalar perturbations sourced by the vector modes 



We are interested in the phenomenological consequences of the vector field production in the model ( 26 ) . As men- 
tioned in the Introduction, we work under the assumption that the 5* — sector has a negligible energy density with 
respect to the inflationary sector, and that the two sectors are coupled only gravitationally. Under these assumptions, 
the main signature of the particle production is encoded in how the gauge quanta enter in the gravitational equations. 
In this subsection we study the effect on the scalar metric and inflaton perturbations; the effect on the tensor modes 
is computed in the next subsection. 

We divide the computation in 3 parts; in the first one, we present the master equation for the density perturbation 
Q; in the second and the third part we compute, respectively, the power spectrum and the bispectrum of the part of 
C, sourced by the vector quanta. 



1. Master equation for 

We perform computations in the spatially flat gauge <5(7ij^scaiar = 0- We then decompose the fields into background 
+ first order + second order perturbations; for example, for the inflaton field we write 

LP = + ^1 + (^2 + ■ • • (37) 

(dots denotes perturbations of higher order, that we ignore) . Under our working assumption that \1/ gives a negligible 
contribution both to the background energy density, and to the cosmological perturbations (which is certainly the 
case, if p^, is sufficiently small), the gauge curvature perturbations C, in this gauge is 

C = -^ ((Pi + (P2) (38) 
9 

We are interested in a master equation for C,. We follow the same steps outlined in Section 5 of |45j. We start from 
the infiaton equation of motion 



5^(^/=55'^''9.^)-\/^^=0 (39) 



Formally, only perturbations of the inflaton and of the metric enter in this equation. However, one can use the 
gravitational equations to express the modes Jgoo and Sg^i in terms of the other perturbations, 

= 5goa [^i, hij] , Sg^i = 6goi A^, /i.y] , (40) 



and then insert these expressions back into (39). The explicit expressions for (40) can be most directly obtained from 
the Einstein equations. ® The vector field modes enter in these equations through their contribution to the energy 
momentum tensor 

T^.,gauge = F^^ " + m2 A^,A, + g^, " ^^') (41) 



where vn? = e^V^°^^ (0- ^ 



We note that only the tensor mode hij enters in the spatial perturbations 5gij, since 5(/ij_scalar = in our gauge, and since we can 
disregard vector metric perturbations as in the standard case. 

We remind that we are working in the unitary gauge, real) and disregarding perturbations of 'I'. In principle, the gauge quanta also 
source perturbations 55', which are then gravitationally coupled to the inflationary sector. We do not expect that this effect is more 
important than the direct gravitational coupling of the gauge modes to the inflationary sector. 
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Equivalently, we can eliminate the Sgoo and Sgoi modes from the action using the so called energy and momentum 
constraints in the ADM formalism. Clearly, this is equivalent to using (40 1 at the level of the equations. Both methods 
are presented in (f 



The net result is that ( 39 ) can be put into a master equation in terms of <^,j 



Ai^„ and h,j. 



The 



master equation, in absence of gauge fields, was first derived in |65j . 

As one can readily employ the same exact steps outlined in [45] to the current model, we do not repeat those 
computations here. There are only two changes that need to be taken into account when one repeats the computations 
of for the present model. Firstly, in the model of the vector field had a direct coupling to the inflaton 
AC — —ffifFF; here the coupling is absent, so we simply formally take j = in all the steps of ^45j. Secondly, the 



4/^ 

vector field was taken to be massless in 



/ 

while it is massive here. This amounts in the additional contributions 
proportional to the vector mass given in (41 ), or equivalently, in the extra term in the cubic action for the perturbations 
in the ADM formalism 



AS': 



1 



(42) 



where A'' is the lapse factor. 



As the energy-momentum tensor (41) is quadratic in the gauge field modes, the first order metric and inflaton 



perturbations are not affected by the gauge field. Therefore the master equation at first order reproduces the standard 
first order expression 



(43) 



where = didi and H = ^. The approximation symbol arises because on the left hand side of this expression we 
have disregarded a "mass term" for ipi which is proportional to the slow roll parameters (the same is true for eq. 
(H) 



M2 



V 



< 1 , 



p Y 



< 1 



(44) 



The reason for this is that we are not interested in slow roll corrections to our results. 
At next order in perturbations theory we find instead 



- V^) + 2H(^^ 



(0)' 



2M2-H 



1 



a''V ■ [Ex B 



(0)' 



2M|,a2 



n 



(45) 



Dots denote terms which are proportional to squares of the first order inflaton and tensor metric perturbations, namely 
they are of O {^p\) and of O (hl ij) (we recall that 6^ can be disregarded); these terms are exphcitly given in [55] . 



The "electric" and "magnetic" fileds are defined in analogy to the electromagnetic expressions, namely 



E, 



1 



A' , B, 



1 



€ijkdjAk 



(46) 



The expression (45) should be compared with eq. (5.15) of [45]. Also here we have disregarded the inflaton "mass 



term", which is given instead in (45j at leading order in slow roll approximation. The expression in [45j had a term 
resulting from the direct inflaton-vector field coupling, that is absent in the present model. On the other hand, 
the final term in ( 45 ) - proportional to the vector mass - is absent in [JS] . Finally, we note the presence of a typo (a 



missing H) in eq. (5.15) of [45j 



2. Scalar source and power spectrum 



Identifying the decompositions ([2]) and (37), the two equations (43) and (45) can then be combined into a unique 
linear non-homogeneous differential equation for Q^, which is formally of the type ([3|, in terms of the two sources 



(47) 



This provides the main difference between the scalar field cosmological perturbations obtained in |45| and in the present work; the 
computations of 1451 were done under the working assumption that the direct interaction is stronger than the gravitational coupling, 
so that all the other terms were disregarded. In this work wo compute instead the more complicated contribution of the gravitational 
interactions. 
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Namely, the source J is obtained from the terms explicitly written on right hand side of (45 1, while J is obtained from 



the terms denoted with dots. The two sources are uncorrelated, and therefore the two particular solutions sourced 
by them can be obtained independently. The solution sourced by J is the one emerging from cosmological second 
order perturbation theory in absence of gauge fields. For standard scalar field slow roll inflation, this term provides a 
negligible contribution to the power spectrum, and unobservable non-Gaussianity. Therefore, we disregard this term 
in this work. In Appendix [B] wc show that 



J: 



2M^Ha 



(27r)^/^ 

where we are using the notation 



(0) 



9(2) fg = Q^f d,g 



(48) 



(49) 



It is worth noting that the source does not diverge as fc — > 0; this is not immediate from (45), since one term contains 
an inverse laplacian. See Appendix [B] for details. 

From (48) we see that, for this model, the source is formally of the type (|4]), with 



r, fc, p] = 



2MjHa 



We insert this operator in (|19|; we also use the fact that I^lpg-j is actually helicity independent in this model, so 



(50) 



that the sum over the helicity a is limited to J2a '^im (p) — ^im —PiPm - see (|14|) - and analogously for the sum over 
a' . Moreover, in the integrand we approximate p— k cri p. We obtain 



fc3 



Gk (r, n) 



Mi a(Ti) 

2 



Gk (r, T2) [ Sp 



dT2 , , 



+ (k-pj 1- (k-py 



l+[k-p 



M2 (ti)M2 (r2)5o"(5g 



k ■ p 



61 SljT^iaM [ri, T2;pr 
(51) 



where we have omitted the sufhx (a) from the correlators (since they are cr— independent in this model). By inserting 
the expressions (34) and (36) for the correlators, and by performing the angular integrals, we obtain 



\\Pp\'[\a,\' + \l3,r 



Gk (t, n) 
a(Ti) 



Gk (r, T2) M(ri)Af (t2) 
a(r2) 



Mrc 



a*pl3p (^i 



,*7(ti) 



27r 

p«7(t2) 



dpp 



11 

15 



Re 



Q,*2^2gi[7(Ti)+7(r2)] _|_ |^^|4gi[7(Ti)-7(T2)] 



(52) 



where 7 (r) = 2 dr' M (r'), and where the Bogolyubov coefficients are evaluated after the particle production, and 



given in ( 29 1 . To obtain this expression we have disregarded p with respect to M inside the mode functions fp present 
inside the correlators. This is a ppropriate since, as we shall see, the integrand is peaked at p '-^ ^/rru, <C M. 

The curly parenthesis in (52) contains several terms proportional to the fast oscillating phase e*^. All these terms 
give a negligible contribution to the final result. To see this, we can separate the momentum and time integrals in 
this expression, and obtain 

PcAk) ^ g^,^2M4 Jdpp'^^\^\'[\ap\' + \(3p\']T^ + ^\l3pnkRe{a;i3p£k)+^^W^^ 



(53) 



The reason for this is that, as we shall see, the signal from particle production is maximal at k ^ H (namely, for modes of the size of 
the horizon when particle production occurs; recall that the scale factor is normalized to one at th e m oment of particle production). 
We shall also see that the source integrand is peaked at p ~ V"^*- We therefore have k <^ p due to teok . 
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where 



dT'Gfc(T, r') m(r') 



T / 

dr Gk [T, T ) m [T ) e ^ > 



Using the expression ( 10 1 for the Green function, the expression 

m (t) = {t — t^) = \n 



H 



— T 



H 



In {-Ht) 



(54) 



(55) 



(in the last equahty we have used the fact that a (r*) = 1) we find, in the super-horizon — fcr ^ 1 hmit (in practice, 
we introduce the integration variable y' — —kr', and we integrate it from to —kr^, rather than from —fcr to —kr^,), 



a (r) TO* 
27 



2F3 



3 3 5 5 5 -fc2\ 
2' 2' 2' 2' 2' 4^2/ 



where the generalized hypergeometric function evaluates to ^ 1 for k <ti H and to - 
The time integral in the exponent of £ can be done analytically, and, leads to 



(56) 

0.423) for k:$> H. 



a (r) 



dy 

-Ht y 



k 

—y cos I —y 



H 



k 

sm I —y 



In (y) I cos 



In y 



\n^y 



(57) 



Since to* ^ H, the term in curly parenthesis is rapidly oscillating unless y = ^ ~ 1. This however means that the 
integrand is peaked at the moment in which the particle production is taking place. In our computation we used the 
asymptotic values (29) for a, /3 (valid once particle production has completed). Therefore we only provide an upper 
bound on Ek- It is easy to verify that instead the integrand of 71- is dominated by times at which (29 1 hold. 

Provided that k <^ \fm^ (which is true in the range of our interest, since the occupation number of the gauge field 
is exponentially suppressed in the opposite regime), we can keep in the integral only the fast oscillating term, and the 
log term, and obtain the estimate 



< 



a (r) TO* 



fc3 
a (r) 
2k^ 



k 
H 



cos 



— sm 



k k k 
_cos--sm- 



for the upper bound. We indeed see that 



< 



dy In (y) 



Ht 



\ny 



I sm 



\ny 



(58) 



< 1. 



This confirms that the oscillatory terms in ( |52[ ) - or (53) 
Performing the momentum integral we finally obtain 



provide a negligible contribution to the final result. 



Pes (fc) 



2 + 72 k^ 



7/2 



46,6567r5 H^M^ 



2F3 



3 3 5 5 5 

2' 2' 2' 2' 2' Alp) 



n 2 



k < 



(59) 



while the result is exponentially suppressed and uninteresting at larger momenta. The exponential suppression of the 
spectrum is due to the fact that the occupation numbers of the vector particles are exponentially suppressed at such 



large momenta, see eq. ( 29 ) . We discuss this result in Subsection III C 



3. Scalar bispectrum 



We now evaluate the formal expression (24) for the bispectrum in this model. This expression contains various 
intermediate quantities that have been given above. Explicitly, we write the green functions in eq. ( 10 ), the projection 
operators in eq. (14), the operator Oip in eq. (50), and the correlators in eqs. (34) and (36). The computation follows 
the same steps presented in the previous Subsection for P^. Also in this case we use the approximation fc ^ p <C M, 
and we find that the terms with a fast oscillating phase in the final time integral can be disregarded. We obtain 



B (fci, k2, fcs) 



27 + 8%/6 



9/2 



22,674,816 7r9/2 H^M^ 



3 3 5 5 5 -fcf 
3 1 2' 2' 2' 2' 2' 4ff2 



(60) 
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This corresponds to the effective equilateral nonlinear parameter (25 1 

9/2 



/jVL,cquil. off. (k) — 2.1 • 10 



2F3 



3 3 5 5 5 -fc2\ 
2' 2' 2' 2' 2' 4ff2 / 



fc < 



(61) 



where we have used the numerical value Pr ~ 2.5 • 10 ^ for the power spectrum. Analogous to the sourced part of the 



power spectrum ( 59 ) , the bispectrum and the nonlinear parameter are exponentially suppressed and uninteresting at 
larger momenta. 



Gravity waves sourced by the vector modes 



We now compute the amount of gravity waves sourced by the produced gauge fields. Inserting the energy momentum 
tensor ( 41 1 for this model in the general expression ^ , we see that the source is formally of the type (|8| with 



n: 



o. 



mn.X 



aJVL 



M' 



^nbjPa {k-p)i 



(62) 



where we recall our notation ( 49 1 



It is instructive to compare this expression with the analogous operator in the source of the scalar perturbations, 
given in ( 50 1 . In writing ( 50 1 we disregarded terms proportional to spatial momenta p and p — k c^i p with respect 
to time derivatives and the mass M. The reason for this is that k p M, as we explained after (52 1. As a 



consequence, time derivatives acting on a mode also give drA ~ MA >• pA, and should be retained in the operator, 
when comp ared to spatial momenta. Disregarding the momentum terms is correct for the scalar source. However, we 
see that in (62 ) the time derivatives enter with an opposite sign to Af^. As a consequence, the dominant contributions 



from the first two terms in (62 1 cancel against each other, and therefore we need to keep the complete structure. We 



now show explicitly how the cancellation arises. 



We insert the operator (62) in the formal expression (21 1 for the power in the tensor modes. The other intermediate 



quantities entering in (21 1 are the Green functions, given in (10), the projection operators, given in (14), and the 
correlators, given in (34) and (36). We obtain 



2fc3 



37r4a2M4 



Gk (t, Ti) 
a(Ti) 



dn 



Gk (t, T2) 
a(T2) 



dpp 



7 r 



Jbl 



^P^ 5bo 



M (ti) SaO - Sal 



- P^ SaO 



2 '4 

M (T2) 5bQ - 6bi + -p SaoS, 



3b0 



T^iafi) [n, T2 \ p] 

(63) 



where we have disregarded k as compared to p. 

This expression is analogous to eq. (51) for Pq, with the difference that here we have already performed the trivial 



angular integrals of d^p. Terms without p in the curly parenthesis are obtained from the square of the first term 



oc 



(^-d^r'' d^r'' +M2^ in (62) ( recall that two Ox enter in Px ) • We note that the structure -di^^d?^ 



is 

,4 



preserved in (63 1 since the suffix (1) on V indicates that A [A') is present in the correlator. The term oc p in 
the curly parenthesis is obtained from square of the other term in (62). The terms (x p^ are the mixed terms. The 
different coefficients (| vs 1) follow from the angular integrals. 



The result (63 1 is a sum of squares of correlators. In each square, most terms present fast oscillating phases; as we 



shall see, these terms give a negligible contribution to the final result once the time integrals are performed. This is 



We remind that the tensor power is obtained from {hxh\) . We have verified that {h\h'^) oc S^x' 
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analogous to |T| ^ \£\ in eq. (54). Squaring the expressions p4| ) and (36 1 we obtain 



2ujp {Tl)ujp (t2) 



+ oscillatory phases 



25(1,0) [n, T2; p]^ = (n) 25(0,0) [n, T2; pf + oscillatory phases 
25(0,1) [ti, T2; p]^ = cjp (T2) 2?(o,o) [n, T2; p]^ + oscillatory phases 
25(1,1) [ti, T2; p]^ cjp (ti) cjp (T2) 2?(o,o) [n, T2; p]^ + oscillatory phases 



(64) 



We recall that ujp (Ti) ~ \J M (r^)^ + p^, and that momentum integrand has its support in the region p <C M {Ti). 

We see, however, that the O (Af^) and O (Af^) parts of the curly parenthesis in (63) cancel for the non-oscillatory 
contributions ( 64 ) . 

Using the full expressions ([34| and ([36]) for the correlators, we obtain 



37r4a2Af4 



2p4 

5wp (Ti)a;p {T2] 



Gk (r, n) 
a(Ti) 

\PpmK\'+\M 



Gk (r, T2) 



a(,T2 
4 



dpp 



p'Wp\ 



(n) 

t^p (Tl)t^p (T2) 



Re 



i7(Ti) 



Tl T2 



2 TM^ (n) (t2) + 6p2 [Af2 (ti) (T2)] + 6p 
5 cjp (ti) (T2) 



Re 



2gi[7(Ti)+7(T2)] 



i[7('ri)-7(i"2)] 



(65) 



where 7 (r^) are the oscillatory phases defined immediately after (52). The first term in the curly parenthesis is 
obtained from the non-oscillatory parts (64), and we see that it is indeed of O (p^). As we shall see, this is the term 
that dominates the final result. 

Eq. ( 65 ) is therefore characterized by a part without oscillatory phase plus a part with one oscillatory phase plus 
a part with two oscillatory phases. For each part we only take the leading prefactor in the p <^ M {Ti) regime. This 
expression then rewrites 



2fc3 



dpp" [\p^\Pp? (|ap|' + \Pp?) fi - Ip'lPpffkRe (a;/3p4) + yRe (^/J^f^ 
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(66) 



where 



Tk 



dY 



, Gk{r,T') 
a (t') ujp (t') 



a{T')ujp (r') 



(67) 



We evaluate the time integrals after the particle production (r > Tmin), so that ap and /3p can indeed be taken as 
constant. Particle production is completed when ujp > ujp, which, in the support region p <, y/rn^ of the momentum 

integral, gives i,nin — -I- yi=- This corresponds to the conformal time Tmin — —jj exp J^- ^ ■ For r > Tmin, we 
can approximate LOp ~ M c^i am, with m given in (55). We then obtain the expression 



Tk 



H^a{T) 
H^a{T) 



cxpl 



Ht 



dy 



H 



H 



cos 



— sm 



y cos 



cxpl 



dy' 
In(y') 



y 



In(y') 
H^a (t) 



sm 



H 



cos 



In 



H 



(68) 



where the integration variable in the first expression is y' = —Ht' . In going from the first to the second line we have 
used the fact that the integrand is peaked at the asymptotically early times, while in the final expression we have 
used H ^ \fml. We note that the result is only logarithmically sensitive to the difference between Tmin and . 
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Under the approximation ut ~ M, the integral £k coincides with the integral £k defined in (54) (we actually notice 
that \£k\ ~ l^fcD- We obtained an upper bound for this quantity in (58). The relative contribution of the two time 
integrals to ( 66 ) can therefore be estimated as 



> In 



IT 



> 1 



(69) 



This shows that the oscillatory integral 8k can be disregarded in our estimate of P^,: 



^ I from At 



157r4a2M4 



Using the results (29), we finally obtain 
P I ~ 



V2) H 



4 . 3/2 

TO* 



32 7r7 M4 fc3 



sm 



k 

cos 

H 



/TO, 

IT 



(70) 



(71) 



while the result is exponentially suppressed and uninteresting at larger momenta. We discuss this result in Subsection 

Unci 

In these last two expressions we have emphasized that in this computation we have considered only the gravity 
waves sourced by the transverse modes of the gauge field. In Appendi x [A| we present the full computation, including 
also the longitudinal vector polarization. We find that the expression ( 70 1 is replaced by 



2P 



Qd At 



157r4a2M4 



/ dpp^ [2|/3,p (|a,p + |/3,n + {\a^\' + W^^] 



(72) 



where and (3^ are the Bogolyubov coefficients of the longitudinal mode. Namely, we see that, if they are produced 
in the same amount, the longitudinal quanta contribute to the gravity waves as one transverse polarization. This 
results in a factor | multiplying the final result (71 ). 



C. Phenomenology 



The observable spectrum of curvature fluctuations in the model ( 26 1 is given by the sum of the sourced contribution 



and the usual nearly scale-invariant contribution from the vacuum fluctuations; see section [Ilj The sourced part of the 
spectrum in the model ( 26 ) leads to a localized "bump" feature in the primordial power spectrum for scales leaving 
the horizon at the moment t — t^,, when the gauge fields are produced. We have illustrated this feature in the left 
panel of Fig. [I] for an arbitrary choice of parameters. The most stringent observational constraints on the model (26) 
come from non-observation of such localized bump features. These observational constraints are the subject of this 
subsection. 

We can write the total observable power spectrum in the following form: 



Pc{k)=V 



(73) 



The spectrum of the vacuum fiuctuations are characterized as usual: kp is the pivot scale (taken to be 0.002 Mpc ^, 
consistent with [66] ) while the amplitude and tilt are given by 



V = 



87r2eM2 



H-2r;-6e. 



(74) 



In this discussion, we disregard the contribution from the longitudinal vector mode to the source of tensor and scalar perturbations. As 
discussed at the end of the previous Subsection, the longitudinal mode changes the result for the gravity wave spectrum by at most a 
factor I (if it is produced in equal amount to each tensor mode). We expect an analogous enhancement for C,. The current discussion 
can be easily modified to account for these additional O (1) factors; all of our conclusions are unchanged. 
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FIG. 1: Left panel: the total observable power spectrum in the model ( 26 \ for a representative choice of parameters, illustrating 



the appearance of a localized bump feature. The solid black curve is the total spectrum while the dashed blue line gives the 
spectrum of the vacuum fluctuations for comparison. Right panel: a comparison of the shape function Si, (solid black curve) 
and the fitting function Sm (dashed red curve). 



Here all quantities are understood to be evaluated at the moment when the pivot scale left the horizon. The sourced 
contribution in (73) describes a localized bump-like feature that we characterize by an amplitude Af,, a location fcf,, 
and a "shape function" Sb- The shape and amplitude are given by: 



Ab = 3.2 • IQ-^V^ 



. 7/2 



Sb[x] = 4.7 c 



2F3 



3 3 5 5 5 
2' 2' 2' 2' 2' 



1 2 



(75) 



The shape function has been constructed so that the global maximum is Sb{x = 1) = 1, meaning that the feature 
in (73 1 reaches its maximum value, Ab, when k = kb- For this work we assume that particle production could have 
taken place at any moment during inflation, so the location of the feature is treated as arbitrary.^'* 

The bump-likc feature ( 75 ) is very similar to the one that would be generated due to instantaneous production 
during inflation of scalar particles directly coupled to the inflaton. The signatures of scalar particle production 
during inflation, including important rescattering effects, were fully derived in [39j . The phenomenology of the bump- 
like feature in the curvature spectrum was subsequently studied in [ID]. (See [H] for a discussion of non-Gaussian 
signatures and [42] for a review.) To derive observational constraints we follow closely the analysis in [40] and replace 



the somewhat complicated shape function ( 75 ) by the following simple fltting function 



which, again, is normalized to have maximal value unity at x — 
formula ( 76 ) provides an adequate description of the feature. 



(76) 



In the right panel of Fig. [T]we show that the simple 



In [40j , a variety of data sets were used to perform a detailed analysis of the observational constraints on localized 
features with shape (76); we refer the reader to that paper for a discussion of the methodology. See also [ST] for 
forecast constraints from Large Scale Structure data, and see ^ST for a discussion of both current and forecasted 
constraints from measurements of the CMB energy spectrum. In the present work we are mostly interested in a 
feature localized on CMB scales, in which case the likelihood contours presented in 00] are approximately flat and 
the observational bound can be roughly summarized as: 



Ab 
V 



1/2 



< 0.1 



H 



< 1.2 7^-1/^ e-2/7 



(77) 



Since the sourced part needs to be subdominant in the power spectrum, V — 2.5 • 10 ^; then, the bound (77 1 can be 
expressed as 



m 



1/2 



H 



< 20e-2/7, 



(78) 



Concretely ft: A.67H if the scale factor is normalized to one at t = 
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1 II 

(Recall that m* ^ H and e <C 1 are both required for theoretical consistency.) We will see shortly that the 
observational bound ( 78 1 excludes the possibility of having any interesting effect from particle production in the 
tensor spectrum. 

So far we have discussed the bump-likc feature in the spectrum of curvature perturbations. However, there will also 
be a corresponding localized feature in the bispectrum. This kind of localized non-Gaussianity is very far from scale 
invariant and hence quite different from the bispectrum templates that are most often used in data analysis. To get 
a very rough sense of the amplitude of non-Gaussianity in our model, we have computed the effective (fc-dependent) 
parameter, which exhibits a bump-like structure. This parameter assumes the maximal value: 



j^equil.cff 
JNL 



(k) 



1.9 • 10" 



-10 



. 9/2 



<92e^/'^ , at A: -1.2/06 



(79) 



It should 
A 



where V — 2.5 • 10 ^ has been used in the first equality and the upper bound comes from imposing (78) 

be emphasized that existing observational constraints (or forecasts) on Z^^"' cannot be applied directly to ( 79 ) 
dedicated search for this type of localized non-Gaussianity would be interesting, however, it is beyond the scope of 
this paper. (See also [HI [55] for a discussion about localized non-Gaussianities from scalar particle production during 
inflation.) In the event of a detection of a bump in the power spectrum, localized non-Gaussian features could play 
an important role to falsify (or support) models of particle production during inflation. 
Using our result (|71|) , the total primordial tensor spectrum in the model ( 26 ) is given by 



where in the shape function 
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0.0226 



sin (4.67 a;) - 4.67 a; cos (4.67 x) 



(80) 

(81) 
0.53, 



the scale x is normalized as in the scalar shape function (75). The shape function is maximized at a; = ^ 
where it evaluates to ~ 1. 

To estimate the amplitude of the gravitational wave signal that can be obtained from particle production, we 
evaluate the tensor-to-scalar ratio on scales where the second term in ( 80 ) is maximized. Assuming the observational 
constraint (78) is satisfied we can approximate ~ Pq^^ = 2.5 • 10" 

3/2 



r{k) 



16e 
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.2 • 10 
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^ and we find: 

1/2^ 



H 



at fc ~ 0.53 kh 



(82) 



Using (78) and e <, 0.006 (corresponding to rvac ^ 0.1) it is straightforward to show that the sourced contribution 



to (82) is always 10 , which is undetectably small. We conclude that an observationally interesting signature in 



gravitational waves cannot be obtained in the model (26). 

Before concluding, we note that the estimates presented here should be interpreted as lower bounds on the efficiency 
of particle production effects. We have considered only a single instance where '0 = 0, leading to a single burst of 
vector particle production. However, in a concrete model one could expect ■0 to undergo damped oscillations about 
the minimum of its potential, passing through zero several times before its kinetic energy is dissipated due to Hubble 
friction or backreaction effects. In such a scenario resonance effects would be expected to enhance the occupation 
number of the produced gauge fields by some factor E that could, in principle, be ^ 1. The enhanceme nt r tfc — >■ En^. 
leads to a factor of E"^ in the sourced contribution to the tensor spectrum, so that the relevant term in (82) becomes 
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(83) 



where we still assume that the scalar spectrum is dominated by th e va cuum fluctuations. The sourced part of the 
scalar spectrum also gets enhanced by a factor of E"^ , so the bound (78) becomes stronger: 

. 1/2 

— < 20(eS)-2/7 . (84) 



H 



Combining these results we find: 



20(£;e 



-2/7 



(85) 



< 10-5 (Etf^ In^ 

max, sourced 

The value E > 200 e"-"^ gives ~ 0.01 at the bump. It would be interesting to study under which conditions this 
value can be reached in a concrete model. 



19 



D. Comparison with GW sourced by modes of different spins 



In Subsection |III B] we have computed the power in gravity waves sourced by vector fields produced in the model 
(26). Schematically, the sourced part of the power spectrum is Vx ex: J d^p (TT), where T is the traceless-transverse 



spatial part of the energy-momentum tensor of source (the gauge fields, in this case) with appropriate contraction. 
The spatial part of the energy momentum tensor contain dominant terms that scale as Tij ^ JVPAiAj in the AI ^ p 
limit (we recall that p and M are, respectively, the momentum and mass of the quanta sourcing the gravity waves). 
One could therefore conclude that V\ oc J dpp^ M*. However, we showed that the dominant terms cancel against each 

other. Also the next to leading term in a Taylor expansion of the integral cancel, and one is left with Vx (x j dpp^, 
see eq. (FtoI). 



In Appendix |D] we performed the analogous computation using fermion fields rather than vector fields as sources. 
In this case the spatial part of the energy-momentum tensor has terms of the type Tij ~ XliPjXy one may 
conclude that V\ oc J dpp'^ M'^ (the factor A/^ coming from the different normalization of the fermion wave function 
with respect to the vector one, compare th e fun ction / in p8| and in (D6l). Also in this case there is however a 



cancellation, resulting in (x^) cx p^, see eq. (D15), and mV\(x J dpp^, see eq. (D16I 



These two scalings agree with that obtained if the source is a scalar particle. In t his c ase, T jj cx piPjc/)^, and one 
immediately has V\ cx J dpp^ without any cancellation. In fact, from our results (70) and (D16), and from the result 
for the analogous computation with a scalar source given in [21], we obtain a very general expression for the power 
spectrum: 



2g. fc3 
157r4a2M4 



f dpp'\P,\'[\ap\' + {-lf' 



(86) 



where s is the spin of the sourcing field, and gs is the number of degrees of freedom of that field: gs = 1 for a scalar, 
gs = 2 for a vector if the longitudinal mode is produced in a negligible amount {g^ = 3 if it is produced in the same 
amount as each transverse mode), and gs = 4: for a Dirac fermion. 

We see that, apart from the difference in the number of degrees of freedom, and a small difference due to the spin 
statistics, the different fields in the nonrelativistic regime M ^ p have a comparable quadrupole moment (transverse 
and traceless projection of Tij) and generate a comparable amount of gravity waves. 



IV. MODEL II: VECTOR PRODUCED BY A PSEUDO-SCALAR INTERACTION 



In this section we consider the following model 



S 



d^x. 



Ml 



--(d^)^~viip)~-id^?-u{ij) 



4/^ 



FF 



(87) 



inflaton sector 



hidden sector 



In addition to a standard inflationary sector, we have introduced a "hidden" sector consisting of a light pseudoscalar, 
-0, and a U{1) gauge field, A^, whose energy density is small a s co mpared to that of the infiaton (so that the Friedmann 



equation takes the usual form 3H « V{(l))). As in section III the hidden sector in (87) has been introduced so that 



the production of gauge field fluctuations can provide a new source of inflationary gravitational waves, complementary 



to the usual quantum vacuum fluctuations of the tensor part of the metric. Unlike the model of section III however, 
we will see that particle production in the theory (87) occurs continuously during inflation, leading to broad-band 



signatures rather than localized features in the scalar and tensor n-point correlation functions. 

The coupling ^p^°\t)FF of the gauge field to the time-dependent pseudoscalar condensate leads to an exponential 
production of fluctuations A^. This effect has already been discussed at length in the literature - see Refs. [i5U45) . 
for example - and here we only review the key features that will be necessary for our analysis. Employing the 
decomposition ([T]) we flnd the following linearized equation of motion of the gauge field mode functions 



2ki 



2Hf 



If the pseudoscalar is in an overdamped regime then the parameter ^ can be treated as a constant. Moreover, we 
assume that tp^^^ > so that the helicity state of the gauge field gets copiously produced while the "— " state 



20 



remains in the vacuum and its effect is renormalized away.^^ The properly normalized solutions of (88) can be written 
as |43l 



1/4 



1/2 



A+{T,k). 



(89) 



This solution is valid only in the phase space interval <C —kr <^ 2^, where the production of gauge fluctuations 
is most important. By restricting ourselves to this regime we effectively cut-off an ultra-violet divergence associated 
with the usual quantum vacuum fluctuations of the gauge field on sub-horizon scales; see [IS] for more discussion. 
We have also assumed that ^ >, 0(1), so that the phase space of produced fluctuations is non-trivial and each mode 
experiences a significant exponential enhancement, e'^^ S> 1, near horizon crossing. (For ^ < 1 there is no interesting 
particle production in the model.) 

We are interested in a scenario where infiation is driven by the potential energy of the (p field, so that 2>H^M^ « Vi^)] 
the "hidden" sector in (87) instead should give a small contribution to the total energy density of the universe. This 



requirement imposes several constraints on the model parameters, which we now discuss. We must first require that 
the energy density in the produced gauge field fiuctuations is smaller than the kinetic energy of 



(90) 



where the "electric" and "magnetic" fields are Ei = —-^A'^, Bi = -^CijidjAi. Using the solution (89) to evaluate the 
expectation value [IS], the condition (90) can be written as: 



(91) 



We also require that the energy density of the rolling pseudoscalar can be neglected with respect to that of the inflaton: 



(92) 



Finally, we require that ^ is adiabatically evolving, = nlpm ~ jp ^ so that it is appropriate to treat it as 
nearly constant during the time interval in which each mode of A is relevant (namely, close to horizon crossing, when 
the mode is produced, and affect cosmological perturbations). As \ ^ 1 during inflation, we need to require 



(93) 



where in the last condition we have approximated U (ip) as a quadratic potential, and we have require d th e evol ution 
of ?/'*-°'' to be in the overdamped regime. Throughout our analysis we will require that the conditions (91), (92), and 
( [93] ) are simultaneously satisfled (these conditions are discussed at the end of Subsection IV C ) . 
The background equation for the pseudoscalar reads: 



+ 3H^(°) + [/' (^(°)) =j{E-B). 



(94) 



It is interesting to note that the condition (91 1 guarantees right hand side of this equation can be disregarded. Indeed 
SSI: 



4^ 



(95) 



which is implied by (91 1 



None of our result for the scalar or tensor correlation functions will depend on the choice i/^^"' > 0. 
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The quantity V^-^^ introduced in ( 13 ) characterizes the two-point function of the produced gauge field fluctuations. 

(96) 



This function, and its derivative, can be written explicitly in terms of the c-number mode functions as 



V 



W 

(0,0) 



[n ,T2;q]=AxiTi, q)Al {t2 , q) , \) [n , T2 ; g] = (n , q)Ai (t2 , q) 



Considering only the "+" helicity state and using the approximate solution ( 89 1 we have 

U/4 



(97) 



A. Scalar perturbations sourced by the vector modes 

1. The master equation 



The hidden sector in (87 1 decouples from the inflaton in the limit Mp — > oo. However, at finite Mp, gravitational 
couplings will transmit the effects of particle production in the hidden sector to the inflaton perturbations, modifying 
the usual predictions for the observable curvature fluctuations. To see this effect we follow closely the analysis of 
subsection III A and derive a master equation for the inflaton perturbations in the model (87). At linear order in 



perturbation theory we have recovered the standard result 



(98) 



where we work to leading order in slow roll parameters. At second order, instead, we have the following equation of 
motion 




(99) 



ivcs J2 



which is formally equivalent to (45) with — >■ 0. In ([99]) the • • • schematically denotes terms involving ip\^ Vi ^^nd 
h\ which do not involve the exponential factors e'^^ that characterize the gauge field modes ( 89 ) and may therefore 
be neglected. 



As we did in Subsection III A both the first and second order equations can be combined into the single master 
equation ([3|. The first order mode ipi is the homogeneous solution of this master equation, while the second order 
(p2 is the particular solution of this master equation. The source in the master equation can be written as a sum of 
two terms 



J^p (t, kj ^ Ji (t, kj + J2 (t, kj 



(100) 



The source J2 is associated with the second term on the right hand side of (99), which is non-local in position space 



We show in Appendix [C| that this term is actually infra-red finite, and that the full source can be cast in the form 

ip~\k-p\r 



J. 





f d^P 


) ~ 4M^H J 


(27r)3/2 



-1 



(r, p) Ei {t, k-pj + Bi (r, p) (T,k-p 



where we have defined the "electric" and "magnetic" field operators as Ei{T,k) 



and Bi 



(101) 



-^eijikjAi yr, kj . This expression appears simpler than the corresponding source (48) in the previous model, due to 

the fact that only the modes A+ are relevant here. 

In this model the energy density in the "electric" field dominates over that in the "magnetic" field [45]. Dropping 
terms involving Bi we have a source term of the form ^ where the model-dependent operator can be approximated 
by 



O. 



{t, k,pj 



'(0) 



ip-\k-p\y 
k^ 



5n di''d?' > 



(102) 



where we recall our notation (49). 
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2. Two-point and three-point correlation functions 



We now proceed to compute the two-point and three-point correlation functions of the gauge invariant curvature 
perturbation, C- We are only interested in contributions that are sourced by particle production effects, since the 
vacuum fluctuations in the model ( 87 1 are standard. The sourced contribution to the power spectrum of ^ is given by 
(19 1. As discussed above, we only consider the ct = 4- contrib utio ns in the sum over helicity states. Using the explicit 

9, a 

dri 



expression (102 1 for the operator O, along with the identity (14), we find that (19) can be written as 
fc3 
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2?;+|)[ti,t2;p] 2?[+i) TuT2;\k~p\ . (103) 



(1,1) 



Next we insert our previous result (97 1 for P^j*^!^. Since we are interested in computing the spectrum at late times, 
the explicit expression ( 10 1 for the Green function can be employed. The sourced power spectrum takes the form 



PcAk) = 



1287r2M4 J (27r) 
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q- jk- q) 
q\k - q\ 



q,\k-q\ 



(104) 



where we have introduced dimensionless variables q = p/k and k = k/k. The dimensionless time integral is defined as 



I [a, b] 



dz 



kr 



smz-zcosz ^2^\V^+Vb] 

^1/2 



(105) 



(notice that z = —kr'). In the super-horizon regime, —kr <^ 1, we can set the lower bound of integratio n to z ero. 
Moreover, particle production effects are most interesting in the regime ^ >, 0(1), in which case the integral (105) has 
most of its support in the region z 1 where we can approximate sin z — z cos z « ^ . Hence, we have the following 
analytical approximation 



I [a, b] 



^5/2 ^ , 

dz e-2y2fi[v^l+Vfc] ^ 
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32\/2 + ^Jb] ^7/2 



(106) 



Finally, the momentum integral in (104) must be performed numerically, giving the final result 

rr4 47r5 
-10 ^ 



PC,.(A:)»4-10- 



M4 e 



(107) 



The computation of the three-point correlation function is completely analogous to the derivation that we have 
outlined for the source power spectrum. Here we simply state the final result for the bispectrum in the equilateral 
configuration: 



Be ^kx =k2^k-i = k)K 2.6 • 10 



_i3 e^'^? 1 



(108) 



The reason for considering the equilateral configuration is that, in this model, the source at any moment is dominated 
by modes with wavelength comparable to the horizon at that moment. This generates mostly correlations between 
scalar perturbations of comparable size [44l |45] . 



B. Gravity waves sourced by the vector modes 



The produced gauge field fiuctuations that are described by the mode solution (89) carry anisotropic stress/energy 
and provide a source of gravitational wave fluctuations that is complimentary to the standard quantum vacuum 
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fluctuations from inflation. -'^^ Tfie computation of the gravitational wave spectrum in the model (87) follows closely 
what we had for the the curvature perturbation. Moreover, this computation has in fact already been performed in 
Therefore we simply state the final result: 



7r2M2 



1 + 8.6-10" 



1 + 1.8 • 10" 



(109) 



The first term in the square braces corresponds to the usual contribution from the quantum vacuum fiuctuations of the 
graviton, and the two helicities have equal power. The second term, on the other hand, corresponds to gravitational 
wave perturbations that have been sourced by particle production effects in the hidden sector. The production is 
much more significant for the /i+ mode, due to the fact that the source consists of only modes. 



C. Phenomenology 



In this subsection we explore the phenomenology of the scalar and tensor cosmological fluctuations in the model 



(87 1. We first consider the spectrum of curvature fluctuations. The total observable power spectrum is the sum of 
(1071 and the standard result for the vacuum fluctuations (see the general result in eq. (19l). Explicitly we have 
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1 + 2.5- IQ-'^e^V- 



V = 



87r2eM2 ' 



(110) 



where the first term in the square braces is the usual spectrum from quantum vacuum fluctuations, while the second 



term is the sourced contribution coming from gauge field production in the model (87). 

Next, we consider the spectrum gravitational wave fluctuations. From (109) and (110) we can write the tensor-to- 
scalar ratio as 



C 



16e- 



l + 3.4-10-5e7'^ 
1 + 2.5 - IQ-^e'^V- 



(111) 



When both the tensor and scalar spectra are dominated by the vacuum fluctuations we recover the standard result, 
r « 16e. On the other hand, if the sourced contributions to (1091 and (110) dominate then we have r k, 218, 
independently of model parameters. The current observational limit is r <, 0.17 [69j . while r ~ 0.01 might be 
detectable with future missions 1 . Since the expression (111) interpolates between 16e and 218, it follows that an 
observable signal can be obtained for any value of e. 



Let us now discuss the phenomenology of the scalar perturbations. From (110) and (111) we get 



Pc,s _ r- 16£ 
P^v ~ 218 



(112) 



We therefore see that the scalar power spectrum is dominated by the vacuum part {Pc^ ~ V). Concerning the scalar 



bispectrum, the effective nonlinearity parameter is immediately obtained from ( 25 ) and ( 108 ) 



-p3 g67r4 
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(113) 



The current CMB bound is 



-214 < < 266 m while /, 



NL 



0(10) might be accessible with Planck [71] 



To obtain the correct amplitude of density fluctuations, we must impose the normalization condition P^ = 2.5 - 10 ^ 
on CMB scales. We can use this condition to eliminate the parameter V in favour of ^ and e. Having done so, the 
key observables r and /^^'' "^ then depend only on the model-dependent quantities ^ and e, which in turn depend 



At second order in cosmological perturbation theory tensor fluctuations can also be sourced by bi-linear combinations of the first 
order scalar fluctuations, 5ip and Sip. Neither of these exhibit the exponential enhancement that characterizes the linear gauge field 
perturbations - see equation (|89| - therefore wc can safely neglect this effect in what follows. 
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FIG. 2: Left panel: The tensor-to-scalar ratio as a function of ^, for several illustrative choices of e. The horizontal line 
corresponds to r = 0.1, the approximate current observational limit. Notice that an observable tensor-to-scalar ratio can 
be achieved for any inflationary potential, by suitably tuning the dynamics in the hidden sector. Right panel: The effective 
nonlinearity parameter as a function of ^, for several illustrative choices of e. The horizontal line corresponds to /jvl ~ 266, 
the approximate current observational limit on non-Gaussianity. 



on the inflationary potential and the dynamics of the hidden sector fields. In Fig. [2] we plot our results for r and /nl 
as a function of ^, for various representative choices of e. 

The observational bound on the tensor-to-scalar ratio forces us into a region of parameter space where non- 
Gaussianity is undetectably small. Therefore gravitational wave fluctuations constitute the most interesting phe- 
nomenology associated with the model (87). This is shown in Fig. [S] where we plot contours in the ^ — e plane leading 
to various phenomenologically interesting scenarios. We note that our findings are relevant also for values of e smaller 
than those shown in the figure. 




FIG. 3: Here we plot contours in the ^ — e plane leading to f^L ~ 266 (the current observational bound on non-Gaussianity), 
r — 0.1 (the current observational bound on tensor modes), and r — 0.01 (which may be detectable in the near future). The 
region above the solid red line is ruled out by producing too much non-Gaussianity while the region above the dashed green 
line is ruled out by non-detection of tensor fluctuations. We see that the non-Gaussianity bound is weaker, meaning that the 
dominant signature of the model comes from gravitational waves. 



As discussed in [33 SSI HI] , the sourced contribution to the tensor spectrum is chiral; only one helicity state is 
efficiently sourced by the gauge field fiuctuations ( 89 ) . This effect may be detected through TB and EB correlations 
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in the CMB [7H [73] . This was first explored by [IB] in the case in which the inflaton is the pseudo-scalar sourcing 
the vector modes; in this case, the direct inflaton-gauge field coupling is so strong that, typically, the main bound 
on the gauge field production is given by the sourced scalar perturbations (non-gaussianity [HI [3S] and, depending 
on the inflaton potential, increased power at small scales [131 [SD]). To overcome this, [IB] assumed the presence 
of ~ 1000 sourcing gauge fields (this decreases the amount of non-gaussianity), or the curvaton mechanism for the 
generation of the scalar perturbations. For some values of parameters, the signal can be above the la detection line 
for a cosmic- variance limited experiment [46j . As we shall now discuss, a more optimistic conclusion is reached if one 
assumes that the gauge field production occurs in a sector only gravitationally coupled to the inflaton, as we have 
studied here. 

A measure of the net handedness of the tensor modes is the following quantity: 
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(114) 



which interpolates between zero (at small ^, when the vacuum fluctuations dominate the tensor mode spectrum) and 
unity (at large ^ when the sourced GW dominate the tensor mode spectrum). In the final approximation we have 
used the fact that, for r < 0.1, the scalar power spectrum in this model is dominated by the vacuum modes. 




FIG. 4: Red/solid lines: Predictions for r vs Ax in the model (871; each line is obtained for a fixed value of e, and for varying ^ 



with greater ^ corresponding to greater particle production, and therefore larger signal. Black/dotted lines: la detection lines 
for the Planck (P), SPIDER (S), CMB-Pol (C), and a cosmic-variance fimited (CV) experiment. The signal needs to be above 
a line to be detectable at la by that experiment. These experimental forecasts are an approximate copy of the lines shown in 
Figure 2 of [73]. 



In Figure|4]we plot the relation ( 114) in the r vs Ax plane, for a few representative values of e; each of the red/solid 



lines is characterized by a given e, and by varying ^ (growing ^ leads to more gravity wave production, and therefore 
greater values of r and £,). We stress that arbitrary large values of r in the range shown in the figure can be reached 
for any value of e. As e decreases, this requires a greater and greater amount of sourced modes, which in turn leads 
to a greater and greater Ax. This explain why, for any given obtained r, greater Ax correspond to smaller e. These 
predictions are superimposed in the figure to la detection lines from various experiments; from top to bottom, the 
lines shown are for the ongoing and forthcoming Planck (P) [TT and SPIDER (S) [71' experiments, for the suggested 
CMB-Pol experiment (C) [T], and for a hypothetical cosmic- variance limited experiment (CV). The signal needs to 
be above a line to be detectable at la by that experiment. These lines arc taken by Figure 2 of [73]. We observe 
that, for some values of parameters, the parity-violation could be detected (at least at la) already by the ongoing / 
forthcoming Planck and SPIDER experiments. 



Before concluding this section, we comment on the constraints (91 1 and ( |92| ) which are necessary for the consistency 
of our calculation. We find: 

0-074^ « « ^ Jl-^ (115) 

^5/2 Mp ^ V V{ip) 
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where the first condition, obtained from (91 1, ensures that the energy density of the produced gauge quanta is smaller 
than the kinetic energy of -0, while the second condition, obtained from (92 1, ensures that the energy density of ^ is 
smaller than that of the inflaton. 

The interval in (115) exists for any reasonable choice of model parameters, therefore we can always choose f/Mp 
such that that various backreaction constraints are satisfied. The lower bound on the inflationary scale becomes 
simpler in the limit in which the sourced part of the gravity wave signal dominates over the vacuum one (which is the 
regime of most interest for our work). Using (111 I in this regime, the lower limit in (115) can be expressed as: 



/ 3.7-10-3^1/4 

w>> — c 



0.01-0.1 



(116) 



where r has been chosen so that the gravity wave signal is observable in the near future (we note that this requires 
C~4-5). 



Finally, using ( 74 ) , the condition ( 93 ) for the adiabatic evolution of ^ can be cast in the form 



Mr 



< 6.6 • 10" 



(117) 



For a quadratic inflaton potential, m^p ~ 6.4 • 10 ^ Mp and ^/e ~ 0.09 (for 60 e-folds of inflation). The condition ( 117) 
then rewrites <^ 9m^. 



V. CONCLUSIONS 



A large experimental effort is currently taking place to detect gravitational waves from inflation. The conventional 
vacuum signal will be detectable only if the scale of inflation is sufRciently high, V^^^ >, 10^^ GeV (corresponding to 
r ^ 0.01, or e >, 6 • 10~^ in single field slow roll inflation). Does this imply that, for a lower scale, this experimental 
effort will be unsuccessful, or can we hope that an observable gravity wave signal can be obtained from some different 
mechanism ? And can we distinguish the gravity waves generated by this mechanism from the conventional vacuum 
ones ? 

Recently, [5T and [25] considered the possible gravity signal from particle production taking place during inflation. A 
difficulty with this idea is that the same particle production will also source scalar perturbations; the observed scalar 
perturbations have properties in perfect agreement with standard vacuum fiuctuations generated by the simplest 
inflationary model; in particular, they are gaussian to a very high degree, while there is no reason to expect that this 
should be the case for a generic model of particle production. Therefore, any mechanism of gravity waves from 
inflation needs to explain why the mechanism has not already manifested itself in the scalar sector. 

To evade this problem, |21| studied scenarios in which the particle production takes place (or increases at a sufficient 
level) only towards the end of inflation, so to generate gravity waves only at scales much smaller than the CMB ones. 
At these scales, the scalar perturbations are only constrained by primordial black holes (see for instance [75] for a 
recent study), which is a significantly weaker bound than those from non-gaussianity at the CMB scales. One of the 
scenarios of [H] was further studied in [33] to account for the backreaction of the produced particles on the background 
inflaton. It was shown in |23| that in this region the sourced gravity wave signal can already be observed by the next 
LIGO stage [4]. The main purpose of the present paper is to study whether an analogous mechanism can produce an 
observable gravity wave signal at CMB scales, without conflicting with the limits from non-gaussianity in the scalar 
sector. 

Let us denote by X the field that is produced during inflation, and that sources the gravity waves. In the models 
considered in [2TJ [22] quanta of X are produced by the motion of the inflaton, which we denote by (p. This implies 
a direct coupling between the inflaton and the produced quanta. As a consequence, if the inflaton is the source of 
cosmological perturbations C, quanta of X will source C with a stronger than gravitational interaction. On the contrary, 
the source of gravity waves from X is of gravitational strength. To minimize the relative amount of produced vs. 
produced gravity waves, in this work we made the opposite assumption of considering the weakest possible coupling 
(in standard gravitational theory) between X and the inflaton: namely, we assumed that ip and X are coupled only 
gravitationally. We therefore assumed that particle production occurs in a "hidden sector" and in this paper, for the 
first time, we computed the amount of scalar perturbations C induced by X through a purely gravitational interaction. 



In fact, as shown in |44| . particle production is a very simple source of observable non-gaussianity of the scalar perturbations, even for 
the highly motivated class of natural inflation models. 
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We then computed the amounts of gravity waves produced by X in these two models (for model II, we quote existing 
results), and we compared the two effects. 

Clearly, there is a large arbitrariness in the choice of the model for particle production, and we do not claim our 
findings to be exhaustive; in particular, we did not study here the analogous of all the scenarios considered in [22) . 
where for instance multiple bursts of particle production, and production of strings were also studied. We study two 
models in which X is a vector field, which is produced by the motion of a field ip ^ (p. In model I, the vector field 
has a mass term -0 (t)^ , and quanta of A are produced when the classical value of ^p crosses zero. In model II the 
vector is continuously sourced by a pseudo-scalar jFF interaction. The reasons for considering these two models 
is that, after the particle production, the vector quanta are highly massive in Model I, while massless in Model II. 
We showed that in the first case this gives rise to a strong suppression of the gravity wave signal with respect to the 
amount of scalar perturbations. In the remainder of this concluding section, we summarize our findings in these two 
models, together with some discussion. 

1. Model I: The main signature of particle production in this model is a bump in the scalar power spectrum, at 
the scales that exited the horizon when the gauge quanta were produced. If this bump will be observed, this 
mechanism can be supported / disproved by the presence/ absence of an analogous bump in the bispectrum, 
which we also computed here for the first time. The spectrum of gravity waves produced by the gauge quanta also 
presents a peak at the same scales, which - if sufficiently high - could distinguish them from the vacuum gravity 
waves. However, we found that - once the bound from not having observed a large bump in the scalar spectrum 
is respected - the amount of gravity waves produced in this model is completely unobservable (r <, 10^®). The 
relative smallness of the gravity waves vs scalar perturbations produced in the model may come as a surprise, 
due to the fact that the gauge quanta are coupled gravitationally to both these quantities. The reason for the 
suppression is due to the fact that the quanta are highly non-relativistic after they are produced. This highly 
suppresses their quadrupole moment, and the amount of gravity waves that they generate. 

To verify this, we also computed the amount of gravity waves produced if the vector field is replaced by a fermion, 
with mass oc "0 (i). Ref. ^1 computed the amount of gravity waves produced by a scalar with mass oc ip {t). From 



our two results, and from the result of [21 for the scalar case, we actually obtained the very general formula (86) 
for the amount of gravity waves produced in all these cases. We see that there is no large enhancement between the 
different spins, apart from the proportionality of the final result to the number of degrees of freedom in each case. 

We did not compute the amount of scalar perturbations C sourced in the fermionic case. However, we believe that 
also in this case the result will be analogous to the one that we have computed, and that therefore our conclusions on 
the relative importance of gravity waves vs scalar perturbations production apply for produced particles of any spin. 

In this model, the vector field is massive after the production, and it therefore also possesses a longitudinal com- 
ponent. The results for this component are more model dependent than those of the transverse components: they 
depend on the specific choice of the potential U for -0 (the results for the transverse component are independent of U 
provided that it is smaller than the kinetic energy of ip during particle production) . The longitudinal mode may drive 
the theory out of perturbative regime when m oc -0 (i) — >■ 0. We showed that this is avoided if the ratio ^ ^ remains 

finite as ^0 — ^ 0- This is, for instance, the case if J7 ~ ^fn^ 0"^ at the origin. These considerations can be relevant 
for all the models in which symmetries are enhanced at some point during the cosmological evolution. For instance, 
we expect massive gauge modes to become massless when different branes move to the same bulk location as in the 
trapping mechanism of [76 . We also found that, for U w ^m'^ip'^, the longitudinal component is produced as much 
as each transverse component, and sources the same amount of gravity waves, for the most reasonable values of the 
mass m^. 

Finally, it is worth pointing out that our study applies to a single instance of particle production. Things may be 
different in cases of multiple bursts of particle production; for instance, if the potential U does not flatten at large ip, 
so that Ip performs oscillations about its minimum, gauge quanta will be produced at each oscillation, in a regime of 
parametric resonance [17j . We have found that the gravity wave signal can reach an interesting level if the parametric 
resonance enhances the amount of produced quanta by a factor of > 200 with respect to the single episode of 
particle production (e being the slow roll parameter). It would be interesting to study under which conditions this 
value can be reached in a concrete model. 



Model II: The amount of gravity waves produced in this model was already computed in [44H46j . The novel 
computation in this work is the amount of scalar perturbations produced in this model under the assumption 



This is qualitatively identical to the findings of | 39| . where the sourced field is a scalar with mass depending on the inflaton. 
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that the inflaton is only gravitationally coupled to the gauge field, and the comparison of the two effects. We 
found that coupling the inflaton only gravitationally sufficiently suppresses the amount of scalar modes generated 
in this model, so that the limits from non-gaussianity are irrelevant when compared to those from gravity waves. 
Therefore, particle production in this model can lead to gravity waves observable at the CMB scales. Differently 
from Model I, the gravity waves in this model are produced at all scales, and not just with a localized bump. 
However, this signal may be distinguishable from the vacuum one since one gravity wave helicity is produced 
in a much stronger amount than the other one, and this can lead |46j to observable nonvanishing TB and EE 
correlations in the CMB 17^175]. 

Ref. |46j already studied whether the parity violation in the sourced gravity waves produced by this model can be 
observed. Also in that case, the problem was to suppress the non-gaussianity of the scalar perturbations produced 
by the gauge field This was overcome in [JS] by assuming the presence of ^ 1000 sourcing gauge fields, or the 
curvaton mechanism for the generation of the scalar perturbations. It was shown in ^46) that, under these assumptions, 
and for some values of parameters, the parity- violation signal can be above the la detection line for a cosmic- variance 
limited experiment }46j. We have seen that in our implementation of the mechanism (namely, by assuming that the 
gauge field is only gravitationally coupled to the inflaton) the parity violation can, for some choice of parameters, be 
observed already by the ongoing / forthcoming Planck and SPIDER experiments. 

We have seen that backreaction bounds from this mechanism are under control for an axion scale / in the interval 
lO^^Mp < / < Mp. Interestingly, the axion decay constant one typically finds in string theory is of the order of 
the GUT scale / ^ 10^^ GeV (see, e.g., [63l [64] 1 which fits comfortably within this window. Indeed, given the UV 
sensitivity of inflation, it is natural to ask whether one can realize our model in string theory. The low energy spectrum 
of string theory contains generically many axion-like particles, which arise from the reduction of antisymmetric p-form 
fields on p-cycles of the internal space. -'^^ These closed string axions have a pseudoscalar coupling i/jFF to U{1) gauge 
fields on the worldvolume of D-branes.^° In such string theory setting, it is not difficult to find inflaton candidates 
with no direct coupling to the axion-gauge field sector. For example, the inflaton can be another axion; the absence 
of direct couplings to the hidden sector follows from some topological and geometrical constraints. Consider a basis 
of p-cycles S^, and their dual p- forms ujj such that ujj = 5j. Two of such axions — J-^ Cp, f = Cp do 

not have kinetic mixing if J uji A *Ldj = 0. The absence of direct inflaton-gauge field ipFF coupling is ensured by 
J uJi Aujj =0 where ujj is the 6—p form dual to ujj, and F is the gauge field to which the axion ip couples. In fact, the 
shift symmetries enjoyed by the axions may provide a natural explanation for why they remain as the light dynamical 



fields during inflation. The modest hierarchy of masses (see e.g., eq. (117)) can be reasonably accommodated without 
necessarily assuming an "axiverse" [75].^^, though having a hierarchy of axion masses offer more flexibilities. Clearly, 
more model building possibilities (with the inflaton being an axion or not) remain to be explored. We hope to return 
to such string theory realizations in the future. 
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Appendix A: Longitudinal vector mode 



In the model (26) the vector has also a longitudinal mode, which was disregarded in the computations performed 



in the main text. The longitudinal sector is actually more subtle than the transverse one, as perturbation theory may 



In addition to these closed string axion-like particles, there are also open string axions but their presence is more model-dependent. For 
this discussion, we shall focus on closed string axions. Moreover, we consider only those axions that are not projected out by discrete 
symmetries (e.g., orientifolding), and do not receive a Stuckelberg mass. 

Such couplings arise from the reduction of Chern-Simons terms in the worldvolume action of D-branes, e.g., ^ Cp A F A F. 

See |80l I81| for some string theory ways of generating a logarithmic hierarchy of axion masses in an "axiverse" . 
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break down in the m — >■ limit. This can be seen, for example, from the action of the longitudinal modes, as we now 
show. 

Let us start by disregarding the perturbations of ^P, as we do in the main text. Also, for simplicity, we disregard the 
expansion of the universe. As we have discussed in the main text, this is a good approximation as long as if ^ \fml. 
Then the mass of the vector field is a classical function of time, and the longitudinal mode is encoded in 



3/2 



A' 



k'lx'l^ - k'Alx' - k^x^Ao + k'\Ao\^ - k^m^lxf + m.'\Ao\' 



(Al) 



The component Aq is non-dynamical (it enters in the action without time derivatives) and it can be integrated out. 
Namely, from (All we obtain the equation 
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2 X 



Plugging this solution back into (|A1|) leads to the action of the longitudinal mode 
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(A3) 



The field L is the canonical field associated with the longitudinal mode. We see that its equation of motion has a 
term that formally diverges as m — 0. The formal reason for the divergence is that the kine tic term for the original 
longitudinal mode vanishes in this limit, as it appears from the relation between L and x (A3 1. This is not surprising, 
since a massless vector has only transverse modes. 

When, as in the present case, the U(l) symmetry is broken spontaneously, there is actually not a decrease of the 
number of degrees of freedom when the classical background part 4'^*^-' vanishes. The physical mass term in the 
original action is obtained from £ D e^|4'(°) (t) + (5^p A^, and the quantity that we have denoted by to in the main 
text is only related to the classical part, rr? = 2e^|4'(°) (i) p. However, when ^^^^ (t) — 0, the fluctuations of ^ 
cannot be disregarded, and one does not obtain a truly massless vector mode at this point (compare this with what 
would happen if to was a hard and time dependent mass in the original theory; in this case there would be instead a 
discontinuity in the number of degrees of freedom at to = 0). 

In short, a full study of the longitudinal sector would require going beyond the linearized theory; this may also affect 
the transverse sector, since all the sectors are coupled to each other beyond the linearized level. Here we use a simpler 
approach, and discuss under which conditions the longitudinal mode does not blow up in the linearized theory; if this 
is the case, then one can expect that dealing with the full theory is unnecessary. We note that these considerations 
apply for a general class of model in which symmetries are enhanced at some point during the cosmological evolution. 
For instance we expect massive gauge modes to become massless when different branes move to the same bulk location 
as in the trapping mechanism of [76]. 



Using the background equations of motion, the dangerous factor m" /m in the linearized computation (A3) rewrites 
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(A4) 



where the first two equations are exact, and in the final step we instead disregard the expansion of the universe. 

Therefore, the equation of motion of L remain finite provided that [/ /ip does not diverge. This could be for 
instance the case for a quadratic potential U ^ |to^V^ for V' ~ 0- In our computations, we assumed a constant 
^(0) ^ rr^ ^ corresponding to V'^"' = ^ (t 
are produced during the time interval At 

potential energy gained by ■f/'^"'' during this interval is smaller than its kinetic energy at (so that ■f/''-''-' can indeed 
be taken as constant) amounts in requiring to^ <C \/m^. This also ensures that the period of oscillations of V'''°' is 



t*). As we discussed in the main text, this implies that the gauge fields 
i — around t — t^. Then, for a quadratic potential, imposing that the 



Clearly, the same equation of motion can also be obtained by writing out the equations for the system in terms of the modes Aq and x 
and by eliminating Aq from these equations. 
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much greater than the time in which particle production takes place, so that one can indeed treat a and /? as constant 
when co mput ing the amount of perturbations sourced by the gauge modes. 

Using (A3), we computed the occupation number of longitudinal vector modes with U ^ ^m'^ip'^, and ^ y/rn^. 
We founoTnat in this regime the longitudinal mode is produced in essentially the same amount as each transverse 
vector mode (for instance, we found that for m^, = 0.1 ^/m.^ the total number densities of the longitudinal and of 
one transverse polarization differ from each other by less than 1%). We stress however that this conclusion is model 
dependent, as the precise evolution of the effective frequency in (A3 1 depends on the details of U (tp). On the contrary. 



the amounts of the transverse modes produced is independent of U, provided that it remains sufficiently smaller than 
the kinetic energy during the time of production. 

In the remainder of this Appendix we study the spectrum of gravity waves produced in this model, once also the 
longitudinal modes are taken into account. When all modes are taken into account, we formally separate the energy 



momentum tensor (41 1 of the vector field into T^^/ + T^J + T^^j", where the first term is quadratic in the transverse 



polarizations, and it is the only one used in the main text, while the third term is quadratic in the longitudinal 
polarization, and the second term is the "mixed term" . Inserting the last two terms into ([t]) we obtain the two 
contributions 
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to the gravity wave source in ([6|. These add up to the term Jj"^, which is the only one studied in Subsection 
of the main text (where it is denoted simply by J\). 

To obtain the total power spectrum we need to evaluate the correlator 

+ + Jx^]ri fci ['^A'^ + ^^A'^ + '^A'T ^2 fc^) and insert it in (11). Each piece in Jx is correlated only 
with the corresponding piece in Jy . 

Proceeding as in the main text, we obtain 
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In these expressions we have introduced the correlators 

(x(ti, Pi)x(t2, P2)> =C(o,o) [n, T2;pi] 6'^^'' {P1+P2) 

which evaluate to 
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where are formally identical to the 'D'^f'J^s quantities given in eqs. ( 34 1 and ( 36 ) , with the only difference that 



in 2?f we use the Bogolyubov coefficients and the mode functions of the longitudinal mode: 



iik) = [aigu+P^gl] af + h.c. 



(AlO) 



31 



where is an annihilation operator, and 
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As in the computation of {J\'^ Jxt^) presented in the main text, the square of the correlators appearing in (A6) 
and ( A7 1 contain terms proportional to fast oscillating phases, which give a negligible contribution to the final result. 
Disregarding these terms, we obtain 
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In these expressions, we have retained only the dominant contributions to the oj^ in the adiabatic regime (namely, we 
have used M' <C Af^ and M" ^ M^). As in the computation of the main text, the terms without oscillatory phases 
are of O (p^) , as a consequence to a cancellation of the would be dominant terms in the transverse and traceless 
projection of the energy-momentum tensor. This cancellation is already visible in (A6| and (A7), see the discussion 
after the analogous expression (63) presented in the main text. 

We see that the mixed term ( A12 1 is suppressed in the adiabatic regime, and can be disregarded. From the (A14) 
we instead obtain 
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where Tk is defined in (67). This contribution adds up to the one of the transverse modes leading to the result (72) 
given in the main text. 



Appendix B: Source for in Model I 



In this appendix, we derive the approximated expression (48 1 for the scalar source in Model I. The source we are 
interested in is the part of (47) that depends on the gauge field. From the three terms at the r.h.s. of (45) we obtain, 
respectively, 
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In the part Ji we have actually disregarded the "magnetic" contribution with respect to the "electric" one. 23 The 

Al 



relative contribution between the two terms in the integrand is -gy ^ ^ ^ 'ap'^ ■ quantity k is the momentum 



An analogous simplification cannot be done for the source of gravity waves, since in this case the dominant term cancels. 
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of the cosmological perturbations that we are computing; we show in the main text that the signal from particle 
production is maximal at k H (we recall that the scale factor is normalized to one at the moment of particle 
production). We also show in the main text that the source integrand is peaked at p ^ ^/m^. We therefore have 
k <^ p due to (30). Finally M ~ am > \fml after the particle production has taken place. As a consequence 

p <^ M and the "magnetic" contribution can indeed be disregarded. 

We perform the time derivative in J2, and eliminate the second derivative through the equation of motion A'^ ~ 
-APA,. We obtain 
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We rewrite ( B4 ) changing the integration variable p 
two. We obtain 



k — p. We add the resulting expression to (B4), and divide by 
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It is worth noting that this expression explicitly shows that J2 does not diverge in the fc — limit. This is not 



immediately obvious from the original expression, since J2 originates from the term with an inverse laplacian in ( 45 1 



Adding this expression for J2 to the expression for Ji and J3 given above, one readily obtains the result (48) given 
in the main text. 



Appendix C: Source for ^ in Model II 



In this appendix, we derive the expression (101) for the scalar source in Model II. In eq. (99) we separated the 
source in the two parts Ji + J2. 



In Fourier space we have the relatively simple expression 

d^p 



Ji 



{r,k) 



AM^nJ (27r)3/2 
for the first part. For the second part, we have instead 
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Some manipulations are instead necessary to put the non-local source term J2 in a form that is amenable to compu- 
tations and is manifestly infra-red finite. We begin by noting that Ai effectively only contains the "+" helicity state 
in this model. This allows us to use the identity 

ki kj ei^^ (p) ef^ (k - pj = p\k - p\ + {k - p) ■ p 

to simplify the tensor structure in the second term of (C2) and re- write J2 as 
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Next we perform the time derivative and use the equation of motion ( 88 ) to eliminate A" which gives the result 
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We rewrite (C5) changing the integration variable p — > k—p. Wc add the resulting expression to (C5), and divide by 
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two. We obtain: 
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We note that this expression is manifestly finite in the limit fc^ — > 0, proving that the non-local source term in (99) 
does not lead to any spurious effects in the infra-red. 

Adding (CI) and (C6I leads to the expression (101) reported in the main text. 
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Appendix D: Fermionic production and gravity waves 

In this appendix, we outline the computation of gravity waves produced by a fermion with a mass varying as in the 
model I studied in the main text. If we denote by X the fermion in the original action, and we rescale, x — X a^^'^, 
the action for the fermion field becomes identical to the one of a massive fermion in Minkowski spacetime, whose mass 
is multiplied by the scale factor 

Sf^ f d^x X [i ~ M (t)] X , M = ga (t) [^(°) (t) - ij;,] (Dl) 



As in the main text, i/j^^^ is a homogeneous classical field that evaluates to at some given moment during 
inflation, while g is the coupling of the Yukawa interaction in the original action. Without loss of generality we 
can choose ijj* — — 0. As in the main text, we consider a regime in which the expansion of the universe can be 
disregarded during the particle production, and we expand 

gil)^°\t)=rn^t , t~0 (D2) 

We instead include the expansion of the universe when we study the amount of gravity waves sourced by the fermionic 
quanta produced at t ~ 0. 

The fermionic production in this model was studied in |77l I78j for reheating after inflation and leptogenesis, and 
in [37] for the imprint on the scalar power spectrum. Here we follow the computations of [78l |82], skipping some 
intermediate steps. We refer the reader to those works for details. One decomposes 
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where -0+ ~ ^ g J ^^'^ ^~ ~ 1 y ' "^^^ spinor V is related to lA by charge conjugation, giving V± (t, fc) = [/^ (r, fc). 

It is convenient to decompose the spinors in terms of the Minkowski solutions (more precisely, the adiabatic solution, 
since in the current case M is not constant) 



t/±(r, A)=a(i)/±(T, fc)T/3W/* (r, fc) , /± (r, fc) ^ ^/l ± ^^e"'^^^^. ''^'^ (D6) 



and E = ^/k'^ + M^. Starting from /3 = at asymptotically early times, and using (D2) one finds, at late times (see, 
for instance, |78) for the computation) 

a(t»0)^^e^^^/V-/V^y^' , /?(t»0)^-e--^V2 ,9^^ (D7) 

As in the bosonic case, is the occupation number, and Pauli blocking is ensured by the fact that, in the fermionic 
case, the Bogolyubov coefficients satisfy |ap + |/?p = l. 

To compute the gravity waves produced by these quanta, we proceed as in the main text, and obtain the formal 
solution (10 1 for the sourced part Qx^s of the canonical gravity wave modes introduced in ([5|. Using the fermionic 
energy momentum tensor 

Tf,^ = ^ [xit^d^X + Xli^df.x - {df,x) l^X - {duX) 7mX] (D8) 



We use 

. -1 

and, only in this Appendix, wc switch to +, — , — , — signature 



— a 
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we can cast the source appearing in ( 10 ) in the form 



1 



2aMr, ^'^'^ 



(2^) 



(D9) 



After Wick contraction, the correlator of two sources acquires the form 

1 



Jx[T,k\jy[T\k 



M^a (r) a (r') 



n 



^ 1 ^hn.X' I ^' 



tr 



(DIO) 



where the transposition acts on the spinor indices. The correlators appearing in this expression need to be regularized. 
We adopt the same prescription adopted in Section [llT] for the vector field correlators. Namely, we normal order the 
fields appearing in the correlators with respect to the time dependent operators 



a{Pj = aa{k^ - f3* (-fc) 



P (-k) = (3a(k 



(Dll) 



that diagonalize the Hamiltonian at any given time. We then recall that the vacuum of the theory is annihilated by 
the original time-independent a and b operators. After some algebra, we obtain 



j/ — p — k 



= ,5(3) (p' ~p^k 
{x (r, P) X (r', k' + ^) y = S(^^ (p- jf - k') 



Cii [r, r'; p'] 1 C12 [t, t'; p'] ? ' P' 
C2i[t, r'; p'] a-p' C22 [t, r'-p'] 1 

Cii + C22 



Cii— C22 , C12+C21 0^ -/ , C12 — C21 ^ ^, 
7 H 77-pH 7-_p 



-C22 [t, t'; p] 1 C21 [r, t';p] a-p 
C12 't'; P] ^-P -Cii [t, t'; p] 1 



(D12) 



where, when particle creation has completed, and (D7| hold, 

Cii [r, r'; p'] = \ [/+ (r,p') /; (r', p') - fl {t,p') /_ (r', p')] - (t,p') /- (r', p') - a*^/! (r,p') /; (r', p')} 

C12 [r, r'; p'] = ^ [/+ (r,p') /* (r', p') + /* (r,p') /+ (r', p')] - (t,p') /+ (r', p') + (r,p') fl (r', p')} 

C21 [r, r'; p'] = -Ci*2 [r, r'; p'] 
C22 [r, r'; p'] = C*i [r, r'; p'] 



(D13) 



We insert these expressions into (DIO). We evaluate the trace, and (as in the analogous computation of the main 
text ) disregard k and fc' in comparison to p and jf . The cPp' integral can be performed using one (5— function. We 
perform the angular part of the remaining cfip integral. Finally, we use the property of the polarization operators 



Ilij. We obtain 



2 5(3) (fc + fc')5;,v 

15 7r2 M^aiT)a{T') 



dpp^ Re [5Cji 



Squaring the correlators, we have 



2E{t) E{t') 



oscillatory phases , = —Cl^ + oscillatory phases 



(D14) 



(D15) 



and we see that the non oscillatory part of the integrand in (D14) is of O (p^) as in the vector case studied in the 
main text. The oscillatory part gives a negligible contribution to the tensor power. Combining this result with (11) 
and ( 20 ) we obtain 



ISTT^a^Mi 



(DI6) 
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where Tk is defined in (67). 

Inserting the result ( |68[ ) for this quantity, and performing the momentum integral, gives 



(8 -~^/2) H'^rr 
167r7M4fc3 



3/2 r 



sm 



k f k 



1 2 



In^ 



k < 



(D17) 



while the result is exponentially suppressed at higher momenta. 
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